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ABSTRACT 


The  subject  of  this  research  is  the  buckling  behavior  of  a  simply  supported 
rectangular  plate,  with  a  bulb-flat  stiffener  attached  to  one  side  of  the  plate.  The  plate 
structure  is  subjected  to  axial  compression  that  increases  to  the  buckling  load.  The 
stiffener  cross-section  has  a  thin  web  and  a  bulb-flat  flange  that  extends  to  one  side  of  the 
web.  Results  of  the  investigation  include  planar  property  formulas  for  the  asymmetric 
flange  geometry,  an  analytic  expression  for  the  Saint  Venant  torsional  constant  of  the 
flange  cross-section,  and  an  analytic  expression  for  the  buckling  load  corresponding  to  a 
tripping  mode  of  the  structure.  The  torsional  constant  for  the  bulb-flat  stiffener  is  15%  - 
23%  higher  than  understood  previously.  The  analytic  expression  for  the  buckling  load  of 
the  bulb-flat  stiffened  plates  considered  in  this  investigation  yields  values  that  are  2%  - 
6%  higher  than  finite  element  results.  It  is  also  shown  that  the  buckling  load  of  a  plate 
with  a  bulb-flat  stiffener  is  3%  -  4%  less  than  that  of  a  plate  with  a  T-flange  stiffener  with 
the  same  cross-sectional  area.  At  the  onset  of  stiffener  tripping,  the  torsionally  superior 
bulb-flat  tends  to  bend  laterally,  while  the  flexurally  superior  T-flange  tends  to  twist. 
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EXECUTIVE  SUMMARY 


Understanding  the  elastic  stability  of  stiffened  plate  structures  is  of  great  interest 
in  the  design  of  bridges,  buildings,  automobiles,  aircraft,  naval  vessels,  and  offshore 
platforms.  This  work  contributes  formulas  intended  for  the  mechanical  engineer,  analyst, 
and  educator  involved  in  the  analysis  or  design  of  the  bulb-flat  cross-section  and  simply 
supported  stiffened  rectangular  plate. 

Chapter  I  comments  on  the  intent  and  synopsis  of  this  investigation.  Brief 
comments  address  the  background  for  the  study  and  literature  search  for  published 
material  related  to  the  investigation.  A  list  of  notation  and  conventions  is  included. 

Chapter  II  presents  the  application  of  multivariable  calculus  to  derive  planar 
property  formulas  for  the  bulb-flat  cross-section.  Appendices  A  and  B  contain  additional 
derivations  and  expressions  related  to  the  formulas.  The  torsional  constant  is  one  of  the 
key  property  values  involved  in  the  analysis  of  stiffened  plate  structures.  No  exact 
expression  exists  in  published  literature  for  the  torsional  constant  of  the  bulb-flat  cross- 
section.  Saint  Venant’s  work  on  an  approximate  expression  for  the  torsional  constant  of 
a  solid  cross-section  provides  a  method  to  develop  an  accurate  one -term  approximate 
formula  applicable  to  bulb-flat  cross-sections  produced  by  Corns  Group,  one  of  the 
world’s  leading  producers  of  bulb-flat  profded  metals.  Presented  is  a  discussion  on  the 
idealization  method  to  approximate  the  torsional  constant. 

Chapter  III  presents  the  application  of  elasticity  theory  and  energy  methods  to 
derive  a  general  expression  for  the  buckling  load  due  to  the  stiffener  tripping  of  a  simply 
supported  rectangular  stiffened  plate  subjected  to  axial  compression.  The  general 
expression  for  the  buckling  load  involves  a  constant  called  ju  .  The  buckling  load  and  ju 
values  are  determined  graphically.  The  value  of  //  indicates  the  stiffener  deflection 
behavior  at  the  onset  of  stiffener  tripping.  Three  special  expressions  provide  initial 
approximations  for  the  upper  bound  of  the  critical  buckling  load  value.  Concepts  applied 
in  the  derivation  include  extremum  principles  in  mechanics,  the  energy  criterion,  calculus 


xv 


of  variations,  the  Rayleigh-Ritz  method,  and  other  topics  in  applied  mathematics  and 
physics.  The  chapter  also  presents  finite  element  analyses  of  several  stiffened  plate 
models.  The  results  are  compared  to  predictions  made  using  derived  fonnulas.  MSC 
Nastran  2001/Patran  2001  r3  finite  element  software  and  the  Maple  8  computer 
environment  are  the  primary  analysis  and  computational  tools  used  in  this  investigation. 
Appendix  C  contains  data  tables  used  to  represent  the  finite  element  models. 

Chapter  IV  summarizes  the  major  findings  and  conclusions  of  this  work.  This 
investigation: 

•  Determines  boundary  equations  for  the  bulb-flat  flange  cross-section. 

•  Derives  planar  property  value  expressions. 

•  Determines  an  approximate  torsional  constant  expression  that  is  more  accurate 
than  idealizing. 

•  Demonstrates  that  the  torsional  property  of  the  bulb-flat  stiffener  is  better  than 
previously  understood.  The  torque-carrying  capacity  of  a  bulb-flat  stiffener 
(possessing  no  structural  flaws)  is  greater  than  that  of  an  area-equivalent  angle 
stiffener. 

•  Derives  a  general  expression  to  predict  the  buckling  load  due  to  stiffener 
tripping  of  a  simply  supported  rectangular  stiffened  plate  subjected  to  axial 
compression.  The  predicted  value  is  less  than  6%  higher  than  the  finite 
element  result. 

•  Demonstrates  that  the  buckling  behavior  of  the  bulb-flat  stiffened  plate  is 
unlike  that  of  the  T-flange  configuration.  The  bulb-flat  tends  to  buckle 
laterally  and  have  a  lower  buckling  load  value  than  an  area-equivalent  T- 
flange  stiffened  plate. 

Chapter  IV  lists  future  research  directions  that  include: 

•  Investigating  methods  to  determine  the  ,u  value  by  other  than  graphical 
means. 

•  Investigating  the  use  of  conformal  mapping  to  determine  the  exact  expression 
for  the  bulb-flat  torsional  constant. 

•  Investigating  solutions  to  the  torsion  problem  for  asymmetric  cross-sections. 


xvi 


Conducting  investigations  of  other  flange  cross-sections,  multiple  stiffener 
configurations,  and  grillages. 

Conducting  further  investigations  of  the  T  flange  by  treating  the  flange  as  a 
thin  web  plate  strip  instead  of  a  beam.  This  treatment  may  achieve  more 
accurate  predictions  for  the  T-flange  buckling  loads. 

Developing  algorithms  that  increase  efficiency  in  stiffened  plate  analysis  and 
design. 
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I.  INTRODUCTION 


A.  INTENT  AND  SYNOPSIS 

This  work  presents  an  investigation  into  the  structural  behavior  and  collapse  of  a 
stiffened  plate  panel  where  the  stiffener  is  thin-webbed  and  the  flange  of  the  stiffener  has 
a  bulb-flat  cross-section.  Chapter  I  includes  brief  comments  on  the  background  and 
motivation  for  the  study.  Chapter  II  presents  analytic  expressions  for  the  cross-sectional 
properties  of  the  bulb-flat  flange.  For  brevity,  Appendices  A  and  B  contain  key 
derivations  related  to  the  property  expressions.  Cross-sectional  properties  include  area, 
centroid,  moments  of  inertia,  and  torsional  constant.  The  property  values  for  the  flange 
cross-section  are  needed  in  the  analytic  fonnulas  to  determine  the  buckling  load  of  the 
stiffened  plate.  Chapter  II  presents  a  comparison  of  the  bulb-flat  stiffener  configuration 
to  the  angle  flange  configuration  to  show  that  the  torsional  property  of  the  bulb-flat 
stiffener  is  better  than  previously  understood.  Chapter  III  presents  a  comparison  of  the 
bulb-flat  stiffened  plate  configuration  to  the  T-flange  configuration  to  show  that  their 
buckling  behaviors  are  different  given  the  same  cross-sectional  area.  These  comparisons 
provide  better  discernment  of  the  buckling  characteristics  of  stiffened  plates.  Chapter  III 
contains  the  development  of  a  buckling  load  formula  for  a  stiffened  plate  resulting  from 
stiffener  tripping.  The  onset  of  stiffener  tripping  negates  the  stiffener’s  support  to  the 
plate  panel  and  leads  to  eventual  collapse  of  the  structure.  The  aim  is  to  develop  an 
analytic  expression  for  the  buckling  load  of  a  rectangular  stiffened  plate,  with  one  or 
more  parallel  bulb-flat  stiffeners  attached  to  one  side,  where  the  plate  structure  is 
subjected  to  axial  compression.  The  focus  of  the  analysis  is  on  a  rectangular  plate 
stiffened  with  one  bulb-flat  stiffener.  The  results  can  be  extended  and  applied  to  plate 
structures  with  more  than  one  stiffener.  Energy  methods  are  used,  as  an  alternative  to 
vector  methods,  to  analyze  the  displacement  and  deflection  behavior  of  the  stiffened  plate 
model.  The  chapter  contains  finite  element  analyses  of  several  stiffened  plate  models 
using  MSC  Nastran  2001/Patran  2001-r3  software.  Finite  element  analysis  results  are 
compared  to  formula  results  to  assess  formula  accuracy  and  the  validation  of  results 
presented  earlier. 
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B.  BACKGROUND 

Stiffened  plates  are  basic  structural  components  of  many  items.  Metal  stiffened 
plates  can  be  found  on  bridges,  buildings,  automobiles,  aircraft,  naval  vessels,  and 
offshore  platforms.  In  ship  design,  a  stiffener  is  a  metal  structure  composed  of  a  web  and 
flange.  Stiffeners  placed  on  one  or  both  sides  of  a  plate  add  strength  and  hinder  overall 
collapse  of  the  plate  panel.  Because  of  the  structural  complexity  of  stiffened  plates, 
understanding  their  elastic  instability  is  of  great  importance  to  ship  designers  and  requires 
careful  study. 

Stiffened  plate  investigations  occur  in  three  general  forms:  theoretical,  using 
classical  and  emerging  theory;  numerical,  using  finite  element  methods  and  computer- 
aided  simulation;  and  experimental,  using  actual  grillages.  It  is  enlightening  to  compare 
the  results  of  all  three  general  forms  of  investigation.  With  regard  to  bulb-flat  stiffened 
plates,  it  appears  that  most  of  what  is  known  of  their  behavior  comes  from  numerical  and 
experimental  investigations.  This  investigation  attempts  to  contribute  new  theoretical 
insights  into  the  behavior  of  bulb-flat  stiffened  plates. 

The  use  of  bulb-flat  plate  stiffeners  in  ship  design  is  said  to  reduce  building  time 
and  maintenance  cost.  Companies  that  produce  bulb-flat  plate  stiffeners  manufacture 
each  stiffener  as  a  single  unit,  which  reportedly  results  in  less  production  cost  compared 
to  welded  or  fabricated  stiffeners.  Additionally,  the  curved  surface  of  the  bulb-flat  traps 
less  moisture  resulting  in  less  corrosion.  The  shape  of  the  bulb-flat  stiffener  is  much 
easier  to  inspect,  weld,  and  paint.  These  benefits  save  significant  repair  and  maintenance 
costs  over  the  lifetime  of  the  vessel. 

Advocates  of  bulb-flat  structures  applaud  the  advantages  gained  by  using  bulb- 
flats.  The  Jiangyin  Bridge  in  China  and  the  Oresund  Bridge  that  links  Denmark  to 
Sweden  are  major  constructions  incorporating  bulb-flat  geometry. 

C.  LITERATURE  SEARCH 

In  his  bibliography,  Langhaar  includes  a  list  of  books  and  articles  that  well  serve 
the  reader  interested  in  general  developments  and  special  topics  in  mechanics.  Notable  is 
the  history  and  theory  of  plate  stability  discussed  thoroughly  in  Bleich  [1],  who 
acknowledges  G.H.  Bryan  as  the  originator  of  the  study  of  plate  stability  under  edge 
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compression.  Bleich’s  discussion  on  rectangular  plates  under  uniform  axial  compression 
and  stiffened  plates  under  axial  compression  are  applicable  to  this  investigation.  His 
demonstration  on  the  theory  of  buckling  of  centrally  loaded  columns  by  torsion  and 
flexure  serves  as  a  fundamental  basis  for  the  buckling  analysis  in  Chapter  III. 

A  series  of  papers  by  Chou  and  Chapman  [3]  are  noteworthy  because  they  address 
the  buckling  behavior  of  bulb-flat  stiffened  plates.  Their  study  presents  an  improved 
design  method,  iterative  in  nature,  for  determining  the  buckling  load  of  the  structure 
against  the  torsional  buckling  of  the  stiffener.  They  validate  the  study  using  FINASIC,  a 
finite  element  program,  and  compare  the  results  to  actual  tests  on  cruciform  struts  and 
box  columns  containing  bulb-flat  stiffeners. 

An  interesting  aspect  of  their  theoretical  analysis  is  the  idealization  of  the  bulb- 
flat  flange  as  an  equivalent  angle  flange.  That  is,  they  treat  the  bulb-flat  cross-section 
like  a  rectangular  cross-section  in  regards  to  the  bulb-flat  torsional  and  warping 
properties.  Chapter  II  addresses  this  treatment  and  shows  that  idealization  imputes  error 
in  the  calculation  of  the  bulb-flat  cross-section’s  torsional  rigidity.  Developing  and  using 
expressions  that  maintain  the  bulb-flat  geometry  yields  more  accurate  cross-sectional 
property  data  resulting  in  a  more  accurate  analysis. 

Other  work  relevant  to  this  investigation  includes  research  and  a  series  of  papers 
by  Danielson  et  al.  [5,  6,  7,  8,  and  9]  analyzing  the  tripping  of  a  beam  attached  to  a  plate 
under  lateral  pressure  loading.  His  investigation  includes  the  assumption  that  the 
stiffener  behaves  like  a  thin-walled  open  section  beam.  Applying  nonlinear  beam  theory, 
he  obtains  an  analytic  expression  for  the  buckling  of  stiffened  plate  structures  under 
longitudinal  compression.  His  work  includes  investigations  of  the  buckling  load  of  ship 
grillages  under  axial  compression  with  and  without  lateral  pressure  [14  and  15].  He  uses 
finite  element  based  eigenvalue  analysis  to  gain  insight  into  the  ways  the  buckling  loads 
and  modes  vary  given  the  grillage  dimensions.  In  new  unpublished  work,  he  extends  his 
previous  work  by  re-deriving  buckling  load  formulas  that  incorporate  different 
assumptions  about  buckling  behavior. 
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D. 


NOTATION  AND  CONVENTIONS 
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Area  of  the  flange 

A, .  +  A  ,  Area  of  the  stiffener 
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Area  of  the  web 
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,  Bending  stiffness  coefficient 


Modulus  of  elasticity  of  the  material  (force/length  ) 
E 


2(l+v) 


Shear  modulus,  modulus  of  rigidity 


Moments  of  inertia  about  the  axes  and  centroid  in  (x,v) 
coordinate  system  (length4) 

Product  of  inertia  about  the  origin  in  (x,  y)  coordinate  system 
(length4) 

Moments  parallel  to  the  axes  through  the  centroid  in  (x,v) 
coordinate  system  (length4) 

Product  of  inertia  about  the  centroid  in  (x,  v)  coordinate  system 
(length4) 

Bending  moments  (force/length)  or  first  moments  (length  ) 

Region  defined  by  the  cross-section  of  the  flange  (length-) 

Height  of  flange  (length) 

Height  of  the  web  (length) 

hw  +  hf ,  Height  of  the  stiffener  (length) 

Slope  of  a  line  or  line  segment 
Radius  of  curvature  (length) 

Thickness  of  the  flange  bulb  (length) 

tw  +  thf ,  Maximum  thickness  of  flange  (length) 
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Thickness  of  plate  (length) 

Thickness  of  the  web  (length) 

Bulb-flat  flange  coordinate  system 

Bulb-flat  flange  coordinate  system  where  the  centroid  coincides 
with  the  origin  of  the  coordinate  system 
Plate  structure  coordinate  system 
Angle  (radians) 

Poisson’s  ratio 

Normal  stress  (force/length“) 
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II.  BULB-FLAT  FLANGE  CROSS-SECTION  PROPERTIES 


A.  INTRODUCTION 

A  wealth  of  knowledge  exists  on  the  behavior  of  stiffeners  with  common 
geometries  such  as  I,  T,  Z,  angle,  channel,  and  flat-bar  cross-sections.  These  common 
geometries  possess  certain  cross-sectional  properties  that  are  relatively  easy  to  derive. 
Studying  uncommon  cross-sections,  such  as  the  bulb-flat  cross-section,  is  often  difficult 
because  obtaining  the  section  properties  of  such  cross-sections  requires  considerably 
more  effort. 

This  chapter  presents  expressions  to  detennine  the  planar  properties  of  the  bulb- 
flat  cross-section.  Subsequent  analysis  of  the  stiffened  plate  structure  requires  knowledge 
of  the  planar  property  values  for  a  specified  cross-section. 

B.  CROSS-SECTIONAL  BOUNDARY  AND  PLANAR  PROPERTIES 
Consider  the  cross-section  of  a  bulb-flat  that  extends  to  one  side  of  the  web  as 

shown  in  Figures  1-3.  Figure  1  shows  the  entire  stiffener  cross-section  composed  of  the 
bulb-flat  flange  and  thin  web.  Figure  2  and  Figure  3  show  only  the  flange  horizontally 
oriented.  Five  independent  variables  ( tw ,  thf ,  r ,  / , ,  and  a)  uniquely  detennine  the 

flange  boundary.  Using  both  Figure  2  and  Figure  3,  the  following  labeling  is  established. 
The  variable  tw  denotes  the  web  thickness  defined  by  the  vertical  distance  between  point 
4  and  point  5.  The  variable  thf  denotes  the  bulb  thickness  defined  by  the  vertical  distance 

between  point  9  and  a  point  horizontal  to  point  5.  The  total  flange  thickness  (at  its 
maximum  value)  is  denoted  by  tf  =  tw  +  thf  defined  by  the  vertical  distance  between  point 

9  and  the  line  formed  by  points  3  and  4.  The  comers  of  the  cross-section  have  radii  of 
curvature  r  and  r\ .  Points  5-6-7  and  8-9-10  define  the  corners  with  radius  of  curvature 

r .  Points  1-2-3  define  the  corner  with  radius  rx .  The  cross-section  has  a  flat  portion 
defined  by  the  points  7  and  8,  that  slopes  a  degrees  from  the  y  -axis. 
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Figure  2  Bulb-flat  flange  geometry 


Figure  3  MSC  Patran  picture  of  bulb-flat  flange  boundary 
points  for  cross-section  properties 
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Let  hf  denote  the  flange  height.  The  boundary  of  the  flange  includes  two  vertical  sides 


x  =  0  and  x  =  h  f ,  where 

hf  =  hf  tana  +  r  (l-  2  tan  a  +  2  sec  a)  (1) 

and  0<a<7T /21  (see  Appendix  A).  The  upper  and  lower  boundaries  of  the  cross- 
section  are  defined  by  the  following  functions.  The  upper  boundary  of  the  flange  cross- 
section  is  given  by: 

fi  (*)  =  (h  +  hf  ~  r)  +  'J1'2  -(x~r)2  ’  0  <x<r  +  r  cos  a 

/(x)  =  <  f2  (x)  =  -x  cot  or  +  tw  +tbf  +  r  (cot  a  +  esc  er-l),  r  +  rcosa<x<hf-rcosa  (2) 

f3  (x)  =  (fi.  +  r)  “  ~(x  ~  hf )  ,  hf  -  r  cos  a  <  x<hf 

The  lower  boundary  of  the  flange  cross-section  is  given  by: 


g(x)  =  <  sAx)  =  h-^rl1-{x-rl)1 , , 
g2(^)  =  °5 


0  <  x  <  i\ 
i\<x<  hf 


(3) 


Given  the  upper,  lower,  and  vertical  boundaries,  the  cross-section  of  the  flange  is 
the  bounded  region  R  defined  by  0  <  x  <  hf ,  g(x)  <  y  <  /(x) ,  with  / (x)  and  g(x)  on 

xe^O,/^.  The  flange  cross-sectional  planar  properties  can  be  calculated  from  the 

following  integral  expression  where  m,n  are  nonnegative  integers  (see  Appendix  B  for 
details  of  the  integral  derivation): 

||  xmy"  dA  =  |  |  x"V  dy  dx  =  — ^  |  xm  {[/(^)]”+1  -  [ g  (^)]"+’}  dx  (4) 

R  0  g(x)  n  0 

Substitution  of  expressions  (2)  and  (3)  into  expression  (4)  yields 


1  As  OC  — >  0 ,  the  flange  geometry  resembles  an  angle  flange  with  a  rounded  end  of  radius  r  whose 
cross-section  properties  could  be  approximated  with  an  equivalent  rectangle. 
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jj  x'”y"dA  = j  (r -l-r  sin#)”'  [tw  +  thf  -r -Hr  cos#)  rcos# 

n  VI  1  77 - 


dO 


hf-rcosa 


where 


I  f  /  \«+i 

+ -  xm  \-xcota  +  tw  +  thf  +  rco\a  +  rQsca-r)  dx 

n  + 1  J  v  w  bf  ' 

r+rcosa 

i  o 

+ - -  J  (hf+r  sin^)  (tw  +  r-r  cos 0)"+'  r  cos 0  dO 

a — 

2 

- l—r”+n+2  J  (l  +  sin^)"'  (l-cos6,),'+1  cos^  dO 


hf  =  thf  tana  +  /'(l-2tana  +  2seca) 

0  <a<~,  0  <x<h 

2  2  2 


/ 


(5) 


Throughout  this  investigation,  the  planar  property  values  for  the  bulb-flat  flange  are 
calculated  using  integral  expression  (5).  Appendix  B  contains  a  summary  of  the  planar 
property  formulas.  For  brevity,  the  expansion  of  each  general  planar  property  formula  is 
omitted. 


Integral  expression  (5)  can  be  specialized  for  a  class  of  bulb-flat  flange  cross- 
sections  by  assigning  a  value  to  one  or  more  of  the  variables.  For  example,  by  setting 
a  =  n  1 6 ,  integral  expression  (5)  defines  planar  property  formulas  for  a  class  of  bulb-flat 
cross-sections  with  a  30°  slope.  Additionally,  if  the  radius  of  curvature  of  the  cross- 
section  comer  rx  is  assumed  to  be  one-tenth  of  the  web  thickness  tw  (/,  =tw/ 10),  the 

integral  defines  planar  property  formulas  for  a  specific  type  of  the  30°  slope  cross- 
sections.  In  this  case,  terms  involving  r\  become  negligible  for  this  particular  class  of 

flange  cross-sections.  Cross-sectional  planar  property  values  for  this  special  class  can  be 
calculated  from  the  following  integral  expression  (6)  where  m,n  are  nonnegative 

integers,  a  =  7l!  6 ,  and  the  tenn  involving  rx  is  neglected: 


11 
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The  following  formulas  are  algebraic  expansions  of  expression  (6)  for  the  specialized 
class  of  bulb-flat  flanges.  The  flange  area  and  centroid  fonnulas  are 
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Formulas  for  the  moments  of  the  specialized  class  of  flanges  are 


l_U/3  2W3  2  13^ 
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4  -  4 + 4  -  4  (^2  +  y2 ) 

C.  TORSIONAL  RIGIDITY 

The  torsion  of  a  beam  is  the  application  of  twisting  moments  at  both  ends  of  the 
beam.  In  the  absence  of  side  surface  tractions  and  if  the  ends  of  the  bar  are  free  to  warp, 
the  beam  experiences  pure  torsion.  The  Saint- Venant  constant  for  uniform  torsion  or 
torsional  constant  J  is  the  ratio  of  an  applied  beam  twisting  moment  Mt  to  the  product 
of  the  shear  modulus  (modulus  of  rigidity)  G  and  the  beam  twist  per  unit  length  0 . 


J  = 


GG 


(7) 


The  role  of  the  quantity  GJ  (torsional  rigidity)  in  the  twisting  of  a  beam  is  similar  to  the 
role  of  the  quantity  El  (flexural  rigidity)  in  the  bending  of  a  beam. 
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1.  Exact  Value  Using  the  Stress  Function  Method 

One  method  of  determining  the  torsional  rigidity  of  a  cross-section  (solving  the 
torsion  problem)  is  by  finding  the  stress  function2.  Analysis  of  torsional  behavior 
requires  finding  the  stress  function  <E> ,  which  is  an  exact  solution  to  Poisson’s  equation: 


32<b  a2o> 


-2 G9  in  R 


dx2  dyJ 

<b  =  0  on  the  boundary  of  R 


(8) 


Donaldson  [10,  p.  376]  shows  the  torsional  constant  as 

R 

The  following  example  illustrates  using  the  stress  function  method  to  determine 
the  torsional  constant.  Let  the  boundary  of  an  elliptic  cross-section  with  the  origin 
coinciding  with  the  centroid  be  given  by 

4  +  £-l  =  0  0°) 

a~  b 

which  has  area  and  inertia  properties 

A  =  nab  Ic=^nab[a2  +Zr) 

Take  the  stress  function  in  the  form 


where  c  is  a  constant.  Hence, 
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(12) 


After  solving  for  c,  equation  (11)  satisfies  the  conditions  of  (8)  and  is  zero  on  the 
boundary  of  the  region  when 


2  Introduced  by  Ludwig  Prandtl. 
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Therefore 
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(13) 


When  a  =  b ,  the  results  correspond  with  the  known  formula  J  =  TTa4  /  2  for  the  circular 
cross-section. 


The  torsion  problem  has  been  solved  exactly  for  other  common  cross-sectional 
forms  such  as  the  square,  rectangle,  and  equilateral  triangle.  The  torsional  constant  can 
be  obtained  relatively  easily  for  flanges  with  cross-sections  that  are  symmetric  about  at 
least  one  axis.  For  flanges  with  arbitrary  cross-sections,  finding  the  exact  solution  is 
quite  difficult.  The  reason  for  this  is,  as  Donaldson  [10,  p.  386]  notes, 

It  is  not  difficult  to  write  expressions  that  satisfy  the  governing  differential 
equation,  and  it  is  not  difficult  to  write  expressions  that  satisfy  the 
boundary  conditions.  What  is  difficult  is  to  do  both  simultaneously. 

The  uncommon  shape  and  nonsymmetrical  cross-section  of  the  bulb-flat  flange  make  it 

difficult  to  find  a  stress  function  that  satisfies  simultaneously  both  conditions  of  (8). 

Several  methods  exist  that  attempt  to  simplify  the  search  for  the  stress  function  of 

“difficult”  cross-sections.3  For  an  arbitrary  cross-section,  the  general  solution  of  the 

torsion  problem  is  found  using  conformal  mapping  by  mapping  the  region  upon  the 

interior  and  boundary  of  a  circle,  then  solving  the  problems  of  Dirichlet  and  Neumann  for 


3  Methods  are  generally  included  in  references  that  discuss  solutions  to  Poisson’s  equation  or  the 
torsion  problem.  Solving  Poisson’s  equation  using  conformal  mapping  is  discussed  in  Henrici  [11,  pp. 
372-377], 
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the  circular  region.  Approximation  methods  offer  alternatives  to  solving  Poisson’s 
equation  directly  for  the  bulb-flat  flange  cross-section. 

2.  Approximate  Expressions:  Lower  and  Upper  Bounds 
Washizu  [24]  uses  variational  methods  and  includes  a  computational  technique 
for  obtaining  accurate  lower  and  upper  bounds  for  torsional  rigidity.  His  method  applies 
the  principles  of  minimum  potential  and  complementary  energy.  Stakgold  [20,  pp.  579- 
583]  presents  a  similar  approach  where  he  derives 


4  JJ  w  dx  dy 


JJ( 


w2  +  w2 )  dx  dy 


-J  -JJ  (Vx2 +Vy2)dxdy 


(14) 


where  w  vanishes  on  the  boundary  of  the  region  and  v  satisfies 

c^v  c^v 

dx2  dy 2 


(15) 


on  the  cross-sectional  region.  The  upper  bound  is  smallest  if  the  region’s  centroid 
coincides  with  the  origin.  Stakgold  further  shows  that  the  bounds  can  be  improved  by 
adding  certain  hannonic  functions  to  v .  It  should  be  noted  that  the  calculations  required 
to  apply  this  method  to  uncommon  cross-sections  is  quite  involved. 

The  follow  examples  illustrate  establishing  bounds  using  (14)  and  (15).  To 
establish  a  lower  bound  for  the  circular  cross-section,  let  w  be  the  following  function 
that  vanishes  on  the  circular  boundary  of  radius  r  whose  center  coincides  with  the  origin. 

w(x,y)  =  x2  +y2  -r2  (16) 
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A  lower  bound  is 


f 

> 

2 

f 

4 

JJ  vv  dx  dy 

4 

JJ  x  +  y2  -r  dx  dy 

R 

\  ^ circle 

7 

V  ^ Circle 

II  ( w  >2  +  w,2 )  dx  dy  JJ  4x’  +  4  v2  dx  dy 


K,rr,r 


For  an  upper  bound,  let 


K,rr.U 


(/  -  Ar 2 ) 


/'  a 

nr 

4 

\2 

-nr 

)  _ 

l  2 

J 

Ur 


Ur 


v(x,y)  =  -Ux2  +  y2) 


(17) 


that  satisfies  condition  ( 1 5).  An  upper  bound  is 

JJ  (v,2  +  v2 )  dx  dy  =  JJ  [x1  +  y2 )  dx  dy  =  / 

circle  < circle 

The  lower  and  upper  bounds  confirm  the  fact  that  for  any  circular  cross-section, 
I  <J  <  I  or  «/  =  /;  that  is,  the  torsional  constant  equals  the  moment  of  inertia  about  the 
centroid. 


For  the  square  cross-section,  with  sides  of  length  2 a ,  establish  a  lower  bound  by 
choosing  w  to  be  the  following  function,  which  vanishes  on  the  boundary  centered  on 
the  origin: 

w(x,y)  =  (x  +  a)(x-a)(y  +  a)(y-a)  =  (x2  ~a2)[y2  -n2)  (18) 

A  lower  bound  is 


u  u 

4  J  J  w  dx  dy 


y 

) 


fl 6a6  V 


IK 


iv,.2  +  wv2  j  dx  dy 


V  9  y 

256 o' 
45 


20  a4 


:  2.2222 a4 


(19) 
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An  upper  bound  is  established  by  letting  v  be  the  same  as  equation  (17). 

a  a  a  a  Q  4 

J  J  (vx2  +  Vv2 )  dx  dy  =  J  J  (x2  +  y2 )  dx  dy  =  — —  =  2.6666a4  (20) 

-a  —a  -a  -a  ^ 

Combining  expressions  (19)  and  (20),  establishes  the  following  lower  and  upper  bounds 
for  the  torsional  constant  of  a  square  cross-section: 

2.2222a4  <J  <2. 6666 a4  (21) 

The  exact  value  of  the  torsional  constant  for  a  square  cross-section  is  J  -  2.2496a4,  as 
shown  later. 

Application  of  the  lower  and  upper  bounds  method  to  the  bulb-flat  flange  cross- 
section  involves  rather  arduous  calculations.  Because  the  upper  bound  is  smallest  if  the 
region’s  centroid  coincides  with  the  origin,  the  boundary  equations  (2)  and  (3)  must  be 
modified  so  that  the  centroid  and  origin  coincide.  Let  the  flange  cross-section  defined  by 
equations  (2)  and  (3)  be  modified  for  the  bounded  region  Rc  (the  subscribe  c  indicating 

the  centroid  coincides  with  the  origin)  defined  by  -x  <  x  <  hf  -  x  ,  g(xc)  <  yc  <  f  (x  ) , 
with  f(xc)  and  g(x)c  on  xc  e  [-x ,  hf  -xj  .  A  function  w  that  vanishes  on  R  is 

wixc’yc)  =  [yc  -  fx  (xc )]  [yc  -f2M][yc-  f3  {*c )] ' x 

(22) 

[yc  -Si  ixc )]  [tc  -  g2  ixc )]  ixc + x)  (xc  -  hf + x ) 

Using  (22)  in  (14)  results  in  very  complicated  expressions  and  provides  only  modest 
gains  to  establishing  accurate  bounds  for  the  bulb-flat  flange  torsional  constant.  Other 
approximation  methods  may  provide  simpler  means  to  obtain  the  torsional  constant  for 
the  bulb-flat  flange  cross-section. 
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3.  Approximate  Expressions:  k*A4/I 

Saint-Venant4  offers  an  approximate  expression  for  the  torsional  constant  of  any 
solid  section  (except  certain  common  sections): 


I  A4 

4x2  Ie 


(23) 


where  A  is  the  cross-sectional  area  and  Ic  is  the  moment  of  inertia  about  the  centroid. 
Let  k  be  defined  as 


kAsL  (24) 

A 

for  a  cross-section  where  expressions  for  the  area,  moment  of  inertia  about  the  centroid, 
and  torsional  constant  are  known.  Consider  the  following  table  showing  the  results  of 
(24)  for  common  cross-sections. 


Table  1  Saint  Venant  approximate  expression  for  common  cross-sections 


Cross- 

Section 

Parameter 

A 

K 

J 

tv 

ii 

Circular 

r 

_  2 
7ir 

—nrA 

2 

1  4 

—nr 

2 

n2r&  1 

,  =  =  0.02533 

4;rV  4  n2 

Elliptic 

a,b 

nab 

^nab^a2  +  b2) 

na2b 3 
a2+b2 

^2  4  r  4  i 

“  4  4  =  2  «  0.02533 

4 nab  4 n 

Square 

2  a 

4  a2 

8  4 

—  a 

3 

2.2496a4 

0.02343 

For  the  given  cross-sections,  the  value  of  k  is  a  fixed  constant  that  does  not  depend  on 
the  parameters  of  the  cross-section. 

In  other  cross-sections,  the  value  of  k  may  vary  as  the  parameters  that  define  the 
cross-section  vary.  For  rectangular  cross-sections  with  length  2b  and  thickness  2 1 ,  k 


4  Discussed  by  Saint-Venant,  “Sur  une  formule  donnant  approximativement  le  moment  de  torsion”, 
Comptes  Rendus,  vol.  88,  1879,  pp.  142-154,  and  in  Timoshenko  and  Goodier  [22,  pp.  301-302]. 
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varies  for  different  values  of  the  ratio  t  /  b  ,5  In  those  cases,  the  torsional  constant  can  be 
obtained  from 


192 1 
7T5b 


t,  4anh 

n= 1,3,5,—  n 


nnb\ 


(25) 


developed  in  Timoshenko  and  Goodier  [22,  pp.  309-313],  For  a  square  cross-section, 
(25)  yields  the  J  expression  shown  for  the  square  cross-section  in  Table  1.  For 
rectangular  cross-sections  with  parameters  2 1  and  2b  ,  the  area  and  inertia  properties  are 


A  =  4bt  Ic=—bt(b2+t 2) 


The  expression  for  k  is  thus 

t=!4=± 

A4  36 


r 


^2  A 


1+- 


V 


V 


1  ,  nnb 

,  .  —  tanh - 

nbn=M,..ns  It 


1  192?  ^ 


Pilkey  and  Chang  [18,  p.  103]  express  the  torsional  constant  fonnula  for  rectangular 
cross-sections  as 


J  =  —  (2h)(2t)3  (1-0.63  — +  0.052— ),  where  t<b  (26) 

which  is  a  simplification  of  (25).  Expression  (26)  gives  adequate  accuracy  and  is  used  in 
the  analysis  presented  in  later  sections.  For  elongated  rectangular  cross-sections,  the 
torsional  constant  can  be  obtained  from 


2  4(2*)(2T 


where  t  <sc  b 


As  the  ratio  t  lb  approaches  zero,  then  (26)  approaches  (27). 


(27) 


5  This  effect  can  be  deduced  from  discussions  found  in  Timoshenko  and  Goodier  [22,  pp.  309-313],  in 
Donaldson  [10,  pp. 390-394],  or  in  Oden  [16,  p.  44]. 
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It  would  be  useful  to  identify  classes  or  sub-classes  of  bulb-flat  cross-sections 
possessing  a  fixed  constant  k  .  Corns  Group6  provides  a  special  profile  brochure  [4]  that 
contains  technical  data  on  59  bulb-flat  cross-sections  commonly  requested.  The 
following  tables  (Table  2  and  Table  3)  reflect  all  59  bulb-flat  cross-sections  published  in 
the  brochure.  The  cross-sections  represent  a  class  of  bulb-flat  flanges  where  a  =  n/6 
and  the  values  for  tw,  tbf,  and  r  are  as  given  in  the  tables.  The  brochure  states  a 

tolerance  level  for  the  value  of  rx .  For  ease  of  calculations,  it  is  assumed  that  rx  =  tw  / 10 , 
which  obeys  the  tolerance  level  stated  in  the  brochure. 


6  Corns  Group  was  formed  in  1999  through  the  merger  of  British  Steel  and  Koninklijke  Hoogovens. 
The  company  is  a  leading  international  metal  company  and  one  of  the  world’s  leading  producers  of  bulb- 
flat  profiled  metals. 
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Table  2  Finite  element  results  for  bulb-flat  flange  cross-sections 


Variables  tw,  thf 

,  and  r, 

are  given,  while  rx 

=  (i/i  OK 

and  a -n ! 6  radians. 

Model 

K 

tbf 

r 

4 

Ic 

J 

k 

mm 

mm 

mm 

mnr 

mm4 

mm4 

1 

6.0 

17.0 

5.0 

335.59 

21903 

13339 

0.023035 

2 

7.0 

17.0 

5.0 

356.15 

24609 

15065 

0.023043 

3 

8.0 

17.0 

5.0 

376.70 

27480 

16915 

0.023084 

4 

6.5 

19.0 

5.5 

410.79 

32907 

19936 

0.023038 

5 

7.0 

19.0 

5.5 

422.19 

34720 

21076 

0.023032 

6 

8.0 

19.0 

5.5 

444.97 

38492 

23475 

0.023049 

7 

10.0 

19.0 

5.5 

490.54 

46667 

28716 

0.023144 

8 

7.0 

22.0 

6.0 

519.50 

53120 

31583 

0.023034 

9 

8.0 

22.0 

6.0 

545.10 

58392 

34823 

0.023031 

10 

9.0 

22.0 

6.0 

570.69 

63911 

38262 

0.023054 

11 

11.5 

22.0 

6.0 

634.66 

78896 

47656 

0.023175 

12 

8.0 

25.0 

7.0 

681.93 

91182 

54641 

0.023039 

13 

9.0 

25.0 

7.0 

711.41 

99088 

59532 

0.023030 

14 

10.0 

25.0 

7.0 

740.89 

107316 

64696 

0.023042 

15 

11.5 

25.0 

7.0 

785.09 

120303 

72912 

0.023088 

16 

8.5 

28.0 

8.0 

849.67 

141249 

85076 

0.023056 

17 

9.0 

28.0 

8.0 

866.36 

146751 

88465 

0.023045 

18 

10.0 

28.0 

8.0 

899.72 

158047 

95505 

0.023035 

19 

11.0 

28.0 

8.0 

933.08 

169751 

102873 

0.023038 

20 

12.0 

28.0 

8.0 

966.43 

181882 

110576 

0.023055 

21 

9.0 

31.0 

9.0 

1035.54 

209455 

126681 

0.023075 

22 

10.0 

31.0 

9.0 

1072.79 

224516 

135976 

0.023049 

23 

11.0 

31.0 

9.0 

1110.03 

240055 

145706 

0.023038 

24 

12.0 

31.0 

9.0 

1147.27 

256095 

155831 

0.023035 

25 

9.5 

34.0 

10.0 

1239.52 

299704 

181896 

0.023095 

26 

10.0 

34.0 

10.0 

1260.02 

309570 

187930 

0.023076 

27 

11.0 

34.0 

10.0 

1301.21 

329720 

200430 

0.023052 

28 

12.0 

34.0 

10.0 

1342.34 

350449 

213450 

0.023039 

29 

10.0 

37.0 

11.0 

1461.61 

416284 

253399 

0.023114 

30 

11.0 

37.0 

11.0 

1506.63 

441894 

269103 

0.023079 
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Table  3  Finite  element  results  for  bulb-flat  flange  cross-sections  (cont) 


Variables  tw,  thf 

,  and  r, 

are  given,  while  /, 

=  (1  / 1 0 )tw  and  a  =  n  /  6  radians. 

Model 

K 

fbf 

r 

4 

K 

J 

k 

31 

12.0 

37.0 

11.0 

1551.64 

468166 

285438 

0.023054 

32 

10.5 

40.0 

12.0 

1701.82 

563855 

344119 

0.023133 

33 

11.0 

40.0 

12.0 

1726.27 

579944 

353915 

0.023113 

34 

12.0 

40.0 

12.0 

1775.17 

612687 

374094 

0.023081 

35 

13.0 

40.0 

12.0 

1824.07 

646209 

395021 

0.023058 

36 

11.0 

43.0 

13.0 

1960.15 

747463 

457230 

0.023151 

37 

12.0 

43.0 

13.0 

2012.93 

787681 

481743 

0.023113 

38 

13.0 

43.0 

13.0 

2065.72 

828771 

507140 

0.023082 

39 

11.5 

46.0 

14.0 

2236.59 

972529 

596102 

0.023167 

40 

12.0 

46.0 

14.0 

2264.93 

997033 

610924 

0.023146 

41 

13.0 

46.0 

14.0 

2321.60 

1046772 

641379 

0.023111 

42 

14.0 

46.0 

14.0 

2378.26 

1097519 

672871 

0.023084 

43 

12.0 

49.0 

15.0 

2531.15 

1244856 

764377 

0.023182 

44 

13.0 

49.0 

15.0 

2591.71 

1304401 

800436 

0.023142 

45 

14.0 

49.0 

15.0 

2652.26 

1365058 

837738 

0.023110 

46 

15.0 

49.0 

15.0 

2712.81 

1426867 

876255 

0.023085 

47 

12.5 

53.5 

16.5 

2990.36 

1736393 

1069001 

0.023213 

48 

13.0 

53.5 

16.5 

3023.55 

1774843 

1092092 

0.023193 

49 

14.0 

53.5 

16.5 

3089.94 

1852703 

1139236 

0.023154 

50 

15.0 

53.5 

16.5 

3156.31 

1931884 

1187887 

0.023123 

51 

16.0 

53.5 

16.5 

3222.69 

2012431 

1237937 

0.023097 

52 

13.0 

58.0 

18.0 

3487.42 

2360220 

1456675 

0.023243 

53 

14.0 

58.0 

18.0 

3559.63 

2458317 

1515332 

0.023202 

54 

15.0 

58.0 

18.0 

3631.84 

2557909 

1575695 

0.023166 

55 

16.0 

58.0 

18.0 

3704.05 

2659049 

1637762 

0.023135 

56 

14.0 

62.5 

19.5 

4061.35 

3199092 

1977224 

0.023249 

57 

15.0 

62.5 

19.5 

4139.39 

3322389 

2050964 

0.023209 

58 

17.0 

62.5 

19.5 

4295.45 

3574247 

2204506 

0.023145 

59 

20.0 

62.5 

19.5 

4529.52 

3966020 

2449749 

0.023082 

] 

Mean  Value 

0.023099 

Standard  Deviation 

5.95E-5 

Median  Value 

0.023084 
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From  Table  2  and  Table  3,  values  for  k  appear  “nearly”  fixed  for  the  specified 
class  of  bulb-flat  flange  cross-sections.  When  applied  to  the  bulb-flat  flange  cross- 
section,  the  Saint  Venant  approximation  (23)  with  k  =  1/  4 n2  ~  0.0253  generates  torsion 
values  that  are  9%  to  10%  larger  than  values  obtained  from  finite  element  methods. 
However,  based  on  the  above  tables,  the  following  approximate  expression  yields 
reasonably  accurate  values  for  the  torsional  constant  of  flanges  within  the  specified  class 
when  compared  to  finite  element  data. 

A  4 

J,  «  0.0231^-  (28) 

7  i, 

Hence,  expression  (28)  serves  as  a  working  fonnula  that  approximates  the  torsional 
constant  for  the  bulb-flange  cross-sections  within  the  specified  class  and  is  used 
throughout  the  remainder  of  this  investigation. 

4.  Idealization 

Success  in  determining  the  exact  expression  of  Saint  Venant’ s  torsional  constant 
is  limited  to  simple  cross-sections.  It  is  common  practice  to  use  known  exact  expressions 
for  simple  cross-sections  to  approximate  more  difficult  and  multiply  connected  cross- 
sections  by  idealization.  Idealizing  is  attributing  a  cross-section,  whose  properties  may 
be  partially  known,  with  the  known  properties  of  a  simple  cross-section.  Caution  is 
required  when  applying  this  method  to  idealize  the  bulb-flat  flange  as  an  area-equivalent 
angle  flange.  The  following  illustrates  mathematically  the  need  for  caution. 

Consider  a  solid  circular  cross-section  of  radius  r  with  area  and  torsional 
properties,  respectively 


A 

^circle 


■  nr  J 


nr 


circle 


An  area-equivalent  square  cross-section  with  sides  2/  in  length  has  area  and  torsional 
properties,  respectively 

4quare=4/2  4quare=2.2496/4 
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Since  AcircIe  =  Asquare ,  then 


/  = 


r 


Hence, 


■Aqua*  =  2.2496/  =  2.2496 


fA  ' 

T" 

_  2.2496;r 

f  nr*  ^ 

2  J 

8 

K  2  y 

0.8834  (yeircle) 


This  shows  that  idealizing  a  circular  cross-section  as  an  area-equivalent  square  cross- 
section  reduces  the  Saint  Venant  torsional  constant  of  the  cross-section  by  nearly  12%. 

A  similar  argument  shows  that  idealizing  an  elliptical  cross-section  as  an  area- 
equivalent  rectangular  cross-section  reduces  the  Saint  Venant  torsional  constant  of  the 
cross-section  when  both  cross-sections  have  high  aspect  ratios.  The  aspect  ratio  for  a 
rectangle  is  the  length  to  thickness  ratio  /  / 1 .  A  high  aspect  ratio  means  the  cross-section 
is  very  long  compared  to  its  thickness.  The  aspect  ratio  for  an  ellipse  depends  on  the 
ratio  of  the  major  axis  length  to  the  minor  axis  length.  Consider  an  elliptical  cross- 
section  centered  on  the  origin  with  semimajor  axis  a  and  semiminor  axis  b  having  area 
and  torsional  properties 


^ellipse  TCak  ^ellipse 


na’b2 
a2 +  b 2 


An  area-equivalent  rectangular  cross-section  of  width  2/  and  thickness  2 1 ,  assuming 
t  <sc  / ,  has  area  and  torsional  properties 


A  =Alt  J 

^rectangle  u  rectangle 


16ft3 

^r 


Figure  4  Idealizing  an  ellipse  as  a  rectangle 


As  commonly  done  in  this  case,  assume  the  two  cross-sections  have  the  same  width 
I  =  a.  Since  Aellipse  =  4ectangle ,  then 

71 

Alt  =  nab  =>  4 at  =  nab  =>  t  =  —  b 

4 


Hence, 


—if  2  .  i2\ 

n  a  +b 

12  [  a2  , 


a2  +b2 ; 


—if  2  .  1 2  N 

7T  a  +b  T 

i  ^  2  J  ellipse 

12  ^  a  J 


The  effect  on  the  torsional  constant  due  to  the  idealization  depends  on  the  relationship  of 
the  semimajor  and  semiminor  axis  lengths.  As  the  ratio  of  the  semiminor  axis  to 
semimajor  axis  approaches  zero,  b  /  a  — »  0  ,  then 


n2  f  a2  +b2' 
12  [  a2  , 


— >  —  =  0.8225 
12 


Even  though  methods  exist  to  minimize  or  eliminate  the  potential  loss  of  torque  capacity 
by  idealization  (detailed  discussion  omitted),  the  argument  is  that  application  of  the 
idealization  method  requires  forethought. 

D.  VALIDATION  USING  FINITE  ELEMENT  MODELS 

Analysis  of  several  finite  element  models  using  MSC  Nastran  2001/Patran  2001- 
r3  software  serves  to  validate  the  property  formulas.  The  property  values  from  formulas 
are  compared  to  values  from  finite  element  models.  In  each  finite  element  bulb-flat 
flange  model,  ten  distinct  points  establish  the  boundary  of  the  bulb-flat  flange  cross- 
section.  Straight  lines  and  3-point  arc  curves  are  defined  from  the  boundary  points.  A 
single  trimmed  surface  is  defined  from  the  curves.  The  arbitrary  shape  option  in  the 
software’s  beam  library  is  used  with  the  maximum  allowable  curvature  error  set  at  0.005. 
Maple  8  computer  environment  and  a  handheld  calculator  are  used  to  calculate  formula 
results.  Table  4  presents  planar  property  values  for  several  bulb-flat  flanges  using 
expression  (28),  fonnula  (76)  in  Appendix  A,  and  formulas  (80)  -  (91)  in  Appendix  B. 
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The  table  presents  a  comparison  of  the  formula  values  to  data  obtained  from  MSC 
Nastran/Patran  finite  element  models  of  the  cross-sections. 
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\Min  Pr 


AREA 

A  =  335.5852 
MOMENTS  OF  INERTIA 

11  =12847.42 

12  =  9055.673 
112  =  -3954.662 

TORSIONAL  CONST. 
ABOUT  CENTROID 
J  =  13339.45 
SHEAR  STIFF.  FACTORS 
K1  =0.7832503 
K2  =  0.7559787 
WARPING  COEFF. 
ABOUT  SHEAR  CENTER 
CWA  =  4531 9.1 7 
CENTROID  RELATIVE 
TO  SHEAR  CENTER 
N1 A  =  0.7564681 
N2A  =  0.6324511 
CENTROID  RELATIVE 
TO  ORIGIN 
HORI.  =  8.1 627645 
VERT.  =  9.3835554 
SHEAR  CENTER 
REL.  TO  ORIGIN 
HORI.  =  7.53031 35 
VERT.  =  8.6270876 
ANGLE  FROM  AXIS  1 
TO  MAX.  PRINC.  AXIS 
ALPHA  =  32.1 9341 
EXTERNAL  PERIMETER 
PERIM  =  76.51 2299 


Figure  5  MSC  Patran  graph  and  data  of  a  horizontal  bulb-flat  flange 


2 


\ 

\ 

\ 

\ 

\Min  Pr 


SECTION  PROPERTIES 
AREA 

A  =  335.586 
MOMENTS  OF  INERTIA 

11  =15800.79 

12  =  8397.743 

II 2  =  -3641 .825 
TORSIONAL  CONST. 
ABOUT  CENTROID 
J  =  13973.02 
SHEAR  STIFF.  FACTORS 
K1  =0.6319036 
K2  =  0.1 45071 6 
WARPING  COEFF. 
ABOUT  SHEAR  CENTER 
CWA  =  449058.7 
CENTROID  RELATIVE 
TO  SHEAR  CENTER 
N1 A  =  7.266891 
N2A  =  1 .493367 
CENTROID  RELATIVE 
TO  ORIGIN 
HORI.  =  7.7303667 
VERT.  =  10.266891 
SHEAR  CENTER 
REL.  TO  ORIGIN 
HORI.  =  6.2369995 
VERT.  =  3. 

ANGLE  FROM  AXIS  1 
TO  MAX.  PRINC.  AXIS 
ALPHA  =  22.26711 
EXTERNAL  PERIMETER 
PERIM  =  87.1 75995 

DIMENSIONS 
W=  20.588 
H  =  23. 

T1  =6. 

T2  =  12.474 


Figure  6  MSC  Patran  graph  and  data  of  an  area-equivalent  angle  flange 


27 


Table  4  Comparison  of  flange  values 


Variables  tw,  thf,  and  r,  are  given,  while  rx  =  (1  /  10)tu  and  a -Tt !  6  radians. 

Percent  change  =  — — — xlOO 

*  -  indicates  an  absolute  percent  change  of  less  than  0.01%. 

Source 

K 

hf 

r 

X 

y 

4 

K 

Percent 

change 

Jf 

Percent 

change 

mm 

mm 

mm 

mm 

mm 

mnT 

mm4 

mm4 

Formula 

6.0 

17.0 

5.0 

8.163 

9.384 

335.59 

21904 

* 

13376 

0.28% 

FE(bulb) 

8.163 

9.384 

335.59 

21903 

* 

13339 

* 

FE(angle) 

7.730 

10.267 

335.59 

24199 

10.48% 

13973 

4.75% 

Formula 

8.0 

17.0 

5.0 

8.397 

10.251 

376.71 

27481 

* 

16928 

0.08% 

FE(bulb) 

8.397 

10.251 

376.70 

27480 

* 

16915 

* 

FE(angle) 

8.010 

11.035 

376.71 

30560 

11.21% 

16026 

-5.26% 

Formula 

7.0 

19.0 

5.5 

9.081 

10.608 

422.19 

34722 

* 

21137 

0.29% 

FE(bulb) 

9.081 

10.607 

422.19 

34720 

* 

21076 

* 

FE(angle) 

8.611 

11.581 

422.19 

38450 

10.74% 

21586 

2.42% 

Formula 

10.0 

19.0 

5.5 

9.408 

11.922 

490.54 

46670 

* 

28661 

-0.19% 

FE(bulb) 

9.408 

11.921 

490.54 

46667 

* 

28716 

* 

FE(angle) 

9.000 

12.754 

490.54 

52004 

11.44% 

26362 

-8.20% 

Formula 

7.0 

22.0 

6.0 

10.036 

11.789 

519.51 

53122 

* 

31674 

0.29% 

FE(bulb) 

10.036 

11.789 

519.50 

53120 

* 

31583 

* 

FE(angle) 

9.486 

12.992 

519.51 

58776 

10.65% 

33449 

5.91% 

Formula 

9.0 

22.0 

6.0 

10.287 

12.642 

570.70 

63914 

* 

38339 

0.20% 

FE(bulb) 

10.287 

12.642 

570.69 

63911 

* 

38262 

* 

FE(angle) 

9.784 

13.735 

570.70 

71244 

11.47% 

36965 

-3.39% 

Formula 

8.0 

25.0 

7.0 

11.574 

13.429 

681.94 

91186 

* 

54785 

0.26% 

FE(bulb) 

11.574 

13.428 

681.93 

91182 

* 

54641 

* 

FE(angle) 

10.939 

14.786 

681.94 

100732 

10.47% 

58328 

6.75% 

Formula 

10.0 

25.0 

7.0 

11.829 

14.281 

740.89 

107321 

* 

64856 

0.25% 

FE(bulb) 

11.829 

14.281 

740.89 

107316 

* 

64696 

* 

FE(angle) 

11.243 

15.528 

740.90 

119359 

11.22% 

63541 

-1.79% 

Examination  of  Table  4  reveals  several  noteworthy  aspects.  The  formula  results 
agree  with  the  finite  element  results  for  the  bulb-flat  flange.  The  two  independent 
sources  agree  almost  exactly  on  the  centroid  locations  (x,y) ,  the  flange  area  values  A  f , 

and  the  values  for  the  polar  moments  /  .  Though  not  independent,  the  fonnula  values  for 
the  torsional  constants  J f  obtained  using  the  approximate  expression  (28)  differs  from 
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the  finite  element  values  by  less  than  one  half  of  a  percent.  Such  agreement  supports  the 
validation  of  the  formula  expressions. 

On  the  other  hand,  note  that  the  finite  element  results  for  the  bulb-flat  flange 
differ  noticeably  with  the  finite  element  results  for  an  area-equivalent  angle  flange.  The 
centroid  coordinates  (x,y)  differ  considerably.  The  moments  /  differ  approximately 

10%.  Only  the  flange  area  results  agree  which  is  by  design.  This  strengthens  the  case 
that  using  an  angle  flange  as  an  area-equivalent  to  the  bulb-flat  flange  produces  error. 

To  determine  the  values  for  the  area  and  torsional  constant  for  the  stiffener  that 
includes  the  web  and  bulb-flat  flange  together,  let  the  stiffener  have  total  height  hs .  The 
stiffener  area  formula  is  given  by 

AJ=A,+Aw=Af+(h,-hr)lw  (29) 

The  membrane  analogy  for  uniform  torsion  is  a  useful  tool  for  visualizing  the 
distribution  of  shearing  stress  in  beam  cross-sections  and  provides  justification  for 
approximating  the  torsional  constant  of  the  stiffener.  If  a  membrane  of  constant  thickness 
were  stretched  over  a  beam  cross-section,  fixed  against  the  boundary  of  the  cross-section, 
and  filled  with  a  gas  exerting  nonnal  pressure  on  the  membrane  surface,  the  volume 
under  the  membrane  is  proportional  to  the  torque-carrying  capacity  of  the  cross-section. 
Additionally,  the  torque-carrying  capacity  of  a  cross-section  composed  of  several 
components  is  greater,  but  not  much  greater,  than  the  sum  of  the  components  treated 
separately.  Hence,  the  total  torsional  constant  of  the  entire  stiffener  is  approximated  by 
the  sum  of  the  component  torsional  constants  (27)  and  (28) 

A  4  1 

J,  >Jr+J„=o.mi-f-+-tJ(h,-hf)  (30) 

In  each  finite  element  stiffener  model,  eleven  distinct  points  establish  the 

boundary  of  the  bulb-flat  stiffener  cross-section  (see  Figure  7  and  Appendix  C).  The 

boundary  points  define  straight  lines  and  3-point  arc  curves.  The  curves  define  a  single 

trimmed  surface.  The  beam  library  provides  the  cross-sectional  properties  for  the  finite 

element  models  using  the  arbitrary  shape  option  with  the  maximum  allowable  curvature 

error  set  at  0.005.  Table  5  compares  the  area  and  torsional  constant  calculations  with 
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finite  element  model  results.  Also,  Table  5  compares  the  area  and  torsional  constant 
calculations  to  published  data  from  the  Corns  technical  brochure  [4], 


Figure  7  MSC  Patran  graphs  of  bulb-flat  stiffener  boundary  points  for  cross- 
section  properties 
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^lax  Pr 


AREA 

A  =  932.0546 
MOMENTS  OF  INERTIA 

11  =1329316. 

12  =  23388.33 
112  =  89132.32 

TORSIONAL  CONST. 
ABOUT  CENTROID 
J  =  20608.63 
SHEAR  STIFF.  FACTORS 
K1  =  0.6633503 
K2  =  0.1 51 8261 
WARPING  COEFF. 
ABOUT  SHEAR  CENTER 
CWA  =  7601410. 
CENTROID  RELATIVE 
TO  SHEAR  CENTER 
N1 A  =  -31 .50906 
N2A  =  2.364666 
CENTROID  RELATIVE 
TO  ORIGIN 
HORI.  =  5.2984028 
VERT.  =  72.076157 
SHEAR  CENTER 
REL.TO  ORIGIN 
HORI.  =  2.9337373 
VERT.  =  103.58521 
ANGLE  FROM  AXIS  1 
TO  MAX.  PRINC.  AXIS 
ALPHA  =  176.1135 
EXTERNAL  PERIMETER 
PERIM  =  275.33551 


Figure  8  MSC  Patran  graph  and  data  of  a  bulb-flat  stiffener  model 
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SECTION  PROPERTIES 
AREA 

A  =  932.0578 
MOMENTS  OF  INERTIA 

11  =1340240. 

12  =  28931 .13 
112  =  101280.4 

TORSIONAL  CONST. 
ABOUT  CENTROID 
J  =  18043.08 
SHEAR  STIFF.  FACTORS 
K1  =0.6610681 
K2  =  0.1 392275 
WARPING  COEFF. 
ABOUT  SHEAR  CENTER 
CWA  =  8301 668. 
CENTROID  RELATIVE 
TO  SHEAR  CENTER 
N1 A  =  -35.1 3593 
N2A  =  2.844582 
CENTROID  RELATIVE 
TO  ORIGIN 
HORI.  =  5.61 64331 
VERT.  =  72.231 934 
SHEAR  CENTER 
REL.  TO  ORIGIN 
HORI.  =  2.771 851 5 
VERT.  =  107 .36786 
ANGLE  FROM  AXIS  1 
TO  MAX.  PRINC.  AXIS 
ALPHA  =  175.6094 
EXTERNAL  PERIMETER 
PERIM  =  286. 


Figure  9  MSC  Patran  graph  and  data  of  an  area-equivalent  angle  stiffener  model 
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Table  5  Comparison  of  stiffener  values 


Variables  hs,  tw,  lhf ,  and  r  are  given.  rx  =  (1  / 1 0)/u,  and  a -Tt !  6  radians. 

J-J(FE) 

Percent  change  = - — y — xlOO .  *  -  absolute  percent  change  of  less  than  0.01%. 

Source 

K 

V 

r 

K 

4 

^c-stifffener 

Percent 

change 

^stiffener 

Percent 

change 

mm 

mm 

mm 

mm 

mm2 

mm4 

mm4 

Formula 

6.0 

17.0 

5.0 

120.0 

932.06 

1352719 

* 

20534 

-0.37% 

FE(bulb) 

932.05 

1352704 

* 

20609 

* 

FE(angle) 

932.06 

1369171 

1.27% 

18043 

-12.45% 

Corns 

931.00 

1353400 

0.05% 

15950 

-22.61% 

Formula 

8.0 

17.0 

5.0 

120.0 

1172.00 

1676051 

* 

33894 

-0.91% 

FE(bulb) 

1172.00 

1676034 

* 

34204 

* 

FE(angle) 

1172.00 

1694072 

1.08% 

30413 

-11.08% 

Corns 

1170.00 

1681000 

0.30% 

27730 

-18.93% 

Formula 

7.0 

19.0 

5.5 

140.0 

1242.45 

2448920 

* 

34535 

-0.46% 

FE(bulb) 

1242.44 

2448895 

* 

34694 

* 

FE(angle) 

1242.44 

2475783 

1.10% 

30482 

-12.14% 

Corns 

1240.00 

2448000 

-0.04% 

27080 

-21.95% 

Formula 

10.0 

19.0 

5.5 

140.0 

1662.34 

3206519 

* 

67721 

-1.25% 

FE(bulb) 

1662.33 

3206493 

* 

68575 

* 

FE(angle) 

1662.34 

3236434 

0.93% 

61791 

-9.89% 

Corns 

1660.00 

3215600 

0.28% 

57520 

-16.12% 

Formula 

7.0 

22.0 

6.0 

160.0 

1460.10 

3788414 

* 

47037 

-0.26% 

FE(bulb) 

1460.09 

3788380 

* 

47160 

* 

FE(angle) 

1460.10 

3831830 

1.15% 

41325 

-12.37% 

Corns 

1460.00 

3788600 

* 

36810 

-21.95% 

Formula 

9.0 

22.0 

6.0 

160.0 

1780.03 

4558106 

* 

70991 

-0.65% 

FE(bulb) 

1780.02 

4558068 

* 

71453 

* 

FE(angle) 

1780.03 

4605079 

1.03% 

63240 

-11.49% 

Corns 

1780.00 

4553200 

-0.11% 

57630 

-19.35% 

Formula 

8.0 

25.0 

7.0 

180.0 

1885.81 

6188305 

* 

70468 

-0.26% 

FE(bulb) 

1885.80 

6188242 

* 

80674 

* 

FE(angle) 

1885.81 

6261495 

1.18% 

70710 

-12.35% 

Corns 

1890.00 

6189000 

0.01% 

63520 

-21.26% 

Formula 

10.0 

25.0 

7.0 

180.0 

2245.73 

7289798 

* 

115017 

-0.57% 

FE(bulb) 

2245.72 

7289731 

* 

115672 

* 

FE(angle) 

2245.73 

7368371 

1.08% 

102123 

-11.71% 

Corns 

2250.00 

7290500 

0.01% 

93280 

-19.36% 
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In  Table  5,  each  stiffener  model  contains  four  rows  of  results.  The  four  rows  are 
referred  to  as  a  result  set  for  clarity.  Each  row  is  data  from  a  different  source.  The  first 
row  in  each  result  set  is  data  from  the  formulas  for  the  bulb-flat  cross-section.  The 
second  row  in  each  set  is  data  from  the  finite  element  analysis  for  the  bulb-flat  cross- 
section.  The  third  row  in  each  set  is  data  from  the  finite  element  analysis  for  the  area- 
equivalent  angle  cross-section.  And  the  fourth  row  in  each  set  is  data  from  the  Corns 
technical  data  brochure. 

Comparing  the  first  two  rows  of  each  stiffener  result  set  compares  the  formula 
results  to  the  finite  element  results.  The  independent  sources  agree  almost  exactly  on  the 
stiffener  area  values  As  and  the  values  for  the  polar  moments  7c_stiffener .  The  formula 

values  for  the  torsional  constants  Jstiffcncr  obtained  using  the  approximate  expression  (30) 

differs  from  the  finite  element  values  by  less  than  two  percent. 

Comparing  the  second  and  third  rows  of  each  stiffener  result  set  compares  the 
finite  element  results  of  the  bulb-flat  to  the  finite  element  results  of  the  area-equivalent 
angle  cross-section.  The  two  sources  agree  on  the  stiffener  area  values  and  the  values  for 
the  polar  moments.  However  there  is  a  marked  difference  for  the  torsional  constant 
values.  Comparing  the  second  and  fourth  rows  of  each  stiffener  result  set  compares  the 
finite  element  results  of  the  bulb-flat  to  the  Corns  technical  data.  The  sources  agree  on 
the  stiffener  area  values  As  and  the  values  for  the  polar  moments.  The  torsional  constant 

values  differ  between  16%  and  23%.  The  reason  for  the  difference  appears  to  be  the  lack 
of  accuracy  when  using  the  idealization  method  to  determine  the  torsional  constant  value. 

In  summary,  this  chapter  presents  expressions  used  to  determine  the  planar 
properties  of  the  bulb-flat  cross-section.  Stiffened  plate  analysis  requires  property  value 
accuracy.  Five  independent  parameters,  web  thickness,  bulb  thickness,  two  radii  of 
curvature,  and  the  slope  angle  uniquely  define  the  bulb-flat  flange  geometry.  The 
essential  planar  property  formula  (5)  takes  the  form  of  an  integral.  All  of  the  property 
values  involve  integral  (5).  The  conventional  methods  of  determining  the  torsional 
constant  are  not  viable  due  to  the  asymmetrical  cross-section  of  the  bulb-flat  geometry 
and  the  complicated  equations  that  define  the  cross-section  boundary.  Expression  (28)  is 
developed  from  Saint  Venant’s  one -term  approximate  expression  and  serves  as  a  working 
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formula  that  yields  very  accurate  results  compared  to  finite  element  methods,  thus 
supporting  its  validation.  Using  the  approximate  expression  provides  values  that  are 
more  accurate  than  those  obtained  by  idealizing  the  bulb-flat  flange  as  an  equivalent 
angle  flange.  The  torsional  constant  for  the  bulb-flat  stiffener  is  15%  -  23%  higher  than 
understood  previously.  In  previous  studies,  the  torsional  nature  of  the  bulb-flat  was 
calculated  using  the  torsional  properties  of  an  area-equivalent  angle  flange  stiffener. 
From  Table  5,  the  torque-carrying  capacity  of  the  bulb-flat  stiffener  is  found  to  be  greater 
than  that  of  an  area-equivalent  angle  stiffener. 
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III.  STIFFENED  PLATE  BUCKLING 


A.  INTRODUCTION 

This  chapter  deals  with  the  buckling  behavior  of  a  stiffened  plate.  One  method  of 
analyzing  the  behavior  of  structural  members  is  by  deriving  the  governing  differential 
equation  and  determining  the  exact  solution  using  an  appropriate  differential  equation 
technique.  When  the  solution  is  difficult  or  impossible  to  determine,  as  is  the  case  with 
stiffened  plate  structures,  approximate  methods  of  analysis  are  used.  The  Rayleigh-Ritz 
method  is  an  effective  alternate  method  that  yields  an  approximate  solution.  The  chapter 
presents  an  analytic  expression  for  the  buckling  load  of  a  rectangular  stiffened  plate,  with 
one  bulb-flat  stiffener  attached  to  one  side,  where  the  plate  structure  is  subjected  to  axial 
compression.  The  results  of  the  single  stiffener  analysis  can  be  extended  and  applied  to 
plate  structures  with  more  than  one  bulb-flat  stiffener.  Additionally,  the  analytic 
expression  can  be  used  to  study  the  behavior  of  stiffened  plates  with  different  flange 
geometries. 

B.  ASSUMPTIONS 


The  fundamental  assumptions  are  taken  from  Danielson  [9],  modified  as 
appropriate,  and  stated  as  follows: 


(0 


Each  plate-stiffener  unit  of  width  b  undergoes  an  identical  deformation. 


(ii) 


(iii) 

(iv) 


The  plate  and  web  obey  the  nonlinear  Von  Karman  plate  equations  (see 
Timoshenko  and  Gere  [21]).  The  flange  obeys  the  nonlinear  beam  equations 
derived  by  Bleich  [1], 

The  plate  and  stiffener  material  is  elastic,  homogeneous,  and  isotropic. 

Every  particle  on  the  bottom  surface  of  a  web  undergoes  the  same 
displacement  as  the  corresponding  particle  on  the  top  surface  of  the  plate,  and 
every  line  of  particles  in  a  web  normal  to  the  plate  surface  remains  normal  to 
the  deformed  plate  at  its  surface.  In  other  words,  the  bases  of  the  stiffeners 
are  clamped  to  the  plate. 
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The  prebuckling  displacements  are  less  than  the  maximum  thickness  of  the 
structure  and  independent  of  the  transverse  coordinate. 

The  incremental  buckling  extensional  strains  at  the  midsurface  of  the  plate 
and  web  are  negligible. 

The  incremental  buckling  displacements  may  be  approximated  by  the 
fundamental  harmonic  in  their  Fourier  expansions. 

The  plate  and  web  are  so  thin  that  their  thicknesses  are  negligible  compared  to 
their  width,  length,  and  wavelength  of  deformation. 

The  stiffener  flange  has  a  solid  cross-section  with  a  length  that  is  much  greater 
than  its  largest  cross-sectional  dimension.  Hence,  the  stiffener  flange  can  be 
treated  as  a  beam. 

The  flange  undergoes  only  lateral  bending  and  torsion,  no  vertical  bending. 

C.  ENERGY  PRINCIPLE 

The  energy  method  for  the  solution  of  elastic  stability  problems  is  based  on  an 
extremum  principle  of  mechanics  that  uses  an  energy  criterion  to  characterize  the 
equilibrium  condition  of  the  elastic  system.  A  more  precise  statement  of  the  energy 
criterion  is  made  in  Danielson  [5]  and  is  recapitulated  in  the  following  discussion. 

Let  Pj  denote  the  potential  energy  of  a  system  in  equilibrium  state  /  .  Let  Pu 

denote  the  potential  energy  in  a  neighborhood  of  equilibrium  state  /  .  Let  P\u\ 
represent  the  increment  in  potential  energy  of  the  elastic  system  in  transition  by 
displacement  field  u  from  state  /  to  a  neighboring  state.  In  addition,  let  P  \ u ]  be 
expandable  into  the  following  component  terms  (functionals) 

P\u]  -  Px  [w]  +  P2  [w]  +  P3  [m]  H — 

where  Tj[w]  refers  to  linear  tenns  (functionals)  with  respect  to  u ,  P2  [«]  refers  to 
quadratic  tenns  (functionals)  with  respect  to  u  ,  and  so  forth.  Hence 

Pu  -  Pj  =  P[u\  =  Px  [i/]  +  P2  [w]  +  P3  [w]  h — 


(v) 

(vi) 

(vii) 

(viii) 

(ix) 

(x) 


36 


If  p[u\  >  0  for  all  non-vanishing  neighborhoods  of  displacement  fields  u  ,  then 

the  fundamental  equilibrium  state  /  is  stable.  In  a  practical  sense,  stability  here  means 
that  a  force,  stress,  shock,  or  disturbance  of  the  system  does  not  cause  an  excessive 
change  or  dramatic  departure  of  the  equilibrium  state  to  a  neighboring  state  or 
configuration.  If  there  exists  a  displacement  field  for  which  P\u\  <  0  ,  the  fundamental 

equilibrium  state  /  is  unstable.  Because  the  elastic  system  is  in  equilibrium,  the  linear 
tenn  of  the  potential  energy  increment  must  vanish.  That  is,  equilibrium  requires 
Px  [i/]  =  0 .  The  condition  necessary  for  stability  is  P2  [w]  >  0  .  The  critical  case  of  neutral 

equilibrium  occurs  when  there  exists  a  displacement  field  ux  such  that  P2  \ux  ]  =  0  and 

P2[u^ux]>  0 .  In  this  case,  the  displacement  field  u]  is  the  buckling  mode,  and  the 

value  of  the  load  that  corresponds  to  this  displacement  field  is  called  the  bifurcation- 
buckling  load.  Stated  mathematically, 

0  <P2[u]  (31) 

is  the  criterion  used  to  determine  the  buckling  mode  and  load  given  the  properly  defined 
total  potential  energy  functional  for  the  stiffened  plate  structure. 

D.  SIMPLY  SUPPORTED  RECTANGULAR  STIFFENED  PLATE 

For  the  stiffened  plate  structure,  the  Cartesian  coordinate  system  (x,  ,x2,x, )  is 
adopted  from  Danielson’s  work  to  aid  result  comparison  and  to  minimize  confusion  with 
the  flange  cross-sectional  coordinate  systems  (x,y)  and  (x  ,  vc ) .  Consider  a  simply 
supported  rectangular  plate  of  length  a  ,  width  b  ,  and  thickness  t  ,  with  a  longitudinal 

stiffener  whose  ends  are  also  simply  supported.  The  stiffener  divides  the  width  of  the 
plate  in  halves  and  is  composed  of  a  thin  web  (treated  as  a  plate  element)  of  height  hw 

and  thickness  tw ,  and  a  bulb-flat  flange  (treated  as  a  beam-column  element).  Suppose  the 
stiffened  plate  structure  is  under  axial  compression  due  to  a  uniform  nonnal  stress  a . 
The  quadratic  terms  of  the  total  potential  energy  functional  for  this  single  stiffened  plate 
unit  at  the  instance  of  buckling  can  be  expressed  as 
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(32) 


P2[W ,V \  -  £/plate  +  t/web  +  f7flange  7^Iate  rweb  Tflange 

W  =  W(x1,x2)  and  V  =  V(xl,x3)  are  the  deflection  fields  (modes)  of  the  plate 
and  stiffener,  respectively.  Subscripts  on  W  and  V  denote  partial  differentiation  with 

d2W 

respect  to  the  coordinate  system,  e.g.  Wl2  = - .  In  expression  (32),  Ui  denotes  the 

oxx  ox2 

internal  potential  energy  as  a  result  of  strains  for  element  i ,  and  Ti  denotes  the  potential 

energy  of  the  external  loads  applied  to  element  i .  According  to  Timoshenko  and  Gere 
[21,  pp.  337,  340,  and  350]  the  internal  potential  energy  expressions  due  to  strains  for  the 
plate  and  web  are  taken  as 

b 

i  2  ci 

Url.t=-Dr\  \{(WU+Wvf  -2(l-v)[wnWa-Wn2  ]U,&2  (33) 

1  _b_  0 

2 

i  K  a 

Umb  =  -  D.  \  j  {(Vt ,  +  V„  f  -  2  (1  -  v)  [V,  iF„  -  V„2  ]}  dXldx,  (34) 

The  potential  energy  expression  for  the  external  loads  applied  to  the  plate  and  web  are 
taken  as 


b 

Opiate  =  T  atp  J  j  W2dxxdx 2  (35) 

2 

i  K  a 

rweb  =  -  (Jt,  j  j  V;dxxdx2  (36) 

^  0  0 

Due  to  assumption  (x),  the  internal  potential  energy  expression  for  the  flange  is 


taken  from  Timoshenko  and  Gere  [21,  p.  25,  eq.  1-52  and  p.  265,  eq.  6-30].  The  value  of 
I xc  is  determined  from  expression  (88). 
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Figure  10  Top  view  of  fiber  curvature  during  buckling 

The  potential  energy  of  the  external  load  for  the  flange  is  the  sum  of  the  products 
of  the  external  forces  and  the  displacements  of  their  points  of  application  in  the  direction 
of  the  forces.  As  the  flange  buckles,  the  stresses  on  the  end  surfaces  may  change  to 
a  +  da .  Considering  small  deformations,  the  assumption  is  that  the  end  conditions  are 
such  that  the  work  done  by  da  may  be  neglected  in  comparison  with  the  work  done  by 
a  .  Hence  the  change  in  potential  energy  for  each  fiber  in  the  cross-section  of  the  flange 
is 

dT^„=(odA)Sc  (38) 

Each  fiber  carries  a  load  adA  and  displaces  longitudinally  a  relative  distance  8C .  The 

displacement  is  due  to  two  reasons:  the  curvature  of  the  fiber  and  the  change  in  the 
longitudinal  stress.  Bleich  explains  and  shows  that  the  potential  energy  due  to  the  change 
in  longitudinal  stress  vanishes.  Hence  the  change  in  potential  energy  for  each  fiber  is  due 
to  fiber  curvature  only.  The  potential  energy  of  the  external  load  for  the  flange  is  derived 
here  and  follows  the  outline  in  Bleich  [1,  pp.  126-127]  and  Chajes  [2,  pp.  204-207]. 
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Let  8C  denote  the  relative  fiber  displacement  due  to  the  curvature  of  the  flange  as 
it  buckles  under  axial  compression.  The  distance  8C  is  equal  to  the  difference  between 
the  arc  length  S  of  the  fiber  due  to  bending  and  the  chord  length  a  of  the  fiber.  Hence 

8 c  - S -a 

The  potential  energy  of  the  external  load  for  the  flange  is  obtained  by  integrating  over  the 
cross-sectional  area  of  the  flange.  Thus 

^flange  =  ojj  S/lx.dx,  (39) 

A 

To  determine  8C,  consider  a  cross-section  of  the  flange  at  distance  x,  along  the 
length  of  the  flange.  The  centroid  of  the  cross-section  coincides  with  the  coordinate 
origin. 


Figure  11  Flange  cross-section  with  centroidal  coordinate  system 

Due  to  buckling,  the  point  with  coordinates  (xc,yc)  in  the  cross-section  of  the  flange  will 
change  coordinates  to  ( xc  +  Ax, ,  yc  +  Ayc ) ,  where  Axc  and  A yc  are  functions  of  x1 .  Let 
the  point  of  rotation  of  the  cross-section  have  coordinates  (£,£)  relative  to  the  centroid 
of  the  cross-section.  Movements  of  the  point  of  rotation  define  the  displacement  of 
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points  in  the  cross-section:  w(jCj)  in  the  xc  direction,  v(x, )  in  the  yc  direction,  and 
/?(xj)  counterclockwise  rotation  about  the  point  of  rotation.  Hence 


Axc  cos  P  -  sin  P  xc  -  £  +  w  xc-E, 

_AycJ  |_sin/?  cos  P  \[_yc -^  +  vj  _yc~C_ 

Since  only  small  deformations  are  considered,  simplification  of  the  right-hand  side  gives 

1  -P  xc-£+w  _  xc-E,  =  w-(yc-C)P~vP 
A Yc_  1  \_yc-C+v_  _yc~C_  v+{xc-t;)p+wp _ 

Using  the  Pythagorean  theorem,  and  considering  <^^X(  =  Axc  and  l0 

dxl  dx j 

small,  yields 

^1  2  2  1  2  2 

<yc=5-a  =  |j  1  +  (axc')  +(Ay/)  -lU^=-f  (a*/)  +(au)  dx\ 

The  potential  energy  of  the  external  load  for  the  flange  is  given  by 

^flange  =|CrjJ'|  (a*/)  +  ( Au/ )  dxpX(dVc  (40) 

Using  the  geometrical  relations  for  the  area,  first  and  second  moments  of  the  flange 
cross-section  (assuming  the  centroid  and  coordinate  origin  coincide), 

JJ  dxcdyc  =  Af 

R 

||  xdxdyc  =  ||  ycdxcdyc  =  0 

R  R 

II  (xc  -  py  +  (yc  -C)2  dxcdyc  =  If 

R 

the  quadratic  terms  of  the  potential  energy  due  to  external  loads  are  given  by 

^flange  =  \ \_Af  (v'2  +  w'2)- 2  A  fit p'  +  2 Af£w  &  + 1 f p'2  dxx  (41) 

^  A 
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Movements  in  the  xc  direction  affected  by  vv  ( x, )  are  assumed  negligible,  by  assumption 
(x),  compared  to  yc  direction  movements  perpendicular  to  the  surface  of  the  web 
yielding 

r„„„  =  \  <t}  ( A,v2  - 2  Afr'P + /,/?'= )  dx,  (42) 

The  point  of  rotation  of  the  cross-section  is  taken  to  be  the  point  centered  on  and  at  the 
top  of  the  web.  This  selection  provides  continuity  of  the  motion  at  the  top  of  the  web  and 
the  flange.  The  point  of  rotation  is  located  half  way  between  points  4  and  5  as  shown  in 
Figure  3  found  on  page  9.  Relative  to  the  orientation  of  the  flange,  is  the  vertical 
distance  from  the  centroid  to  the  point  of  rotation.  For  the  bulb-flat  cross-section 

%  =  hf -x  =-(x  -  hf^  (43) 

and 


where  x  is  determined  from  (83)  and y  is  determined  from  (84).  A  lateral  movement  of 
the  point  of  rotation  v(x, )  equates  to  V  (xt,hw).  A  rotation  or  twist  about  the  point  of 
rotation  of  the  flange  /3  [xx )  equates  to  V3  (;q ,  hw ) . 
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1.  Energy  Functional 

The  quadratic  terms  of  the  total  potential  energy  functional  of  a  single  stiffened 
plate  unit  at  the  instance  of  buckling  can  be  expressed  as 


p  =  Tj  +U  +U  -T  -T  -T 

x  2  plate  web  flange  plate  ±  web  flange 


b 

2  a 


11^22  ^12  J  j  dX\dx2 


Z  _b  0 
2 

+  ^DX!{(Vn+v»f-Hl-v)[W»-^]}&A,  +  ^S(EI^+GJ/Va2)  ^  dx,(4S) 


1 

- G 

2 


o  o 

b 

2  a 


J  J  tIflt2dxtdx1  +  J  J  tHP12dxldx3  +  J  [drV]1  +  2 Af  (x-h^  VxVn  +  IfVu2  j  dxx 

_b  0  0  0  0 

2 

The  boundary  conditions  of  the  stiffened  plate  structure  are 


W(0,x2)  =  W(a,x2)  =  0 

(46) 

Wn(0,x2)  =  Wu(a,x2)  =  0 

(47) 

V[0,x3 )  =  V  ( a,x3 )  =  0 

(48) 

Vn  (0,x3)  =  Vn  (a,x3)  =  0 

(49) 

V(Xl,0)  =  W(xx,0) 

(50) 

V3(xl,0)  =  W2(xl,0) 

(51) 

2.  Approximate  Solution 

The  Rayleigh-Ritz  method  is  a  very  effective  technique  that  finds  an  approximate 
solution  to  problems  requiring  the  first  variation  of  the  total  potential  energy  functional  to 
vanish.  The  method  is  summarized  in  two  steps.  First,  assume  an  admissible  solution 
containing  unknown  coefficients  that  satisfy  the  boundary  conditions  of  the  problem. 
Second,  substitute  the  assumed  solution  into  the  functional  and  determine  the  value  of  the 
unknown  coefficients  that  minimizes  the  functional.  From  experimental  observation  and 
by  assumption  (vii),  the  deflection  of  the  plate  panel  is  assumed 
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(52) 


\  .  mnx,  .  7tx1 

xl,x2)  =  csin - Lsin — - 

a  b 

where  c  is  an  arbitrary  constant.  The  value  m  is  the  number  of  half  waves  in  the 
buckling  mode  along  the  longitudinal  axis  and  is  taken  to  be  the  integer  that  gives  the 
lowest  value  for  the  buckling  load.  The  deflection  of  the  stiffener  is  assumed 

V  (x1 ,  x3 )  =  c  sin  m7tXx  /  (x3 )  (53) 

where  / (x,)  is  a  function  to  be  determined.  Expressions  (52)  and  (53)  satisfy  boundary 

conditions  (46)  -  (49).  Boundary  conditions  (50)  and  (51)  transform  to  the  following 
conditions: 


/(o)=o 

(54) 

T(o)=f 

(55) 

Hence, 


1  2  2  _4  /  < 

1  c  m  77  /  1 


2  4  ah1 

Term,  =  Dh^t  f< 

7T  J 


/I 

f  mb  a  ^ 

( ~D,P 

- 1 - 

\2  p 

a  mb  y 

+  Tenn,  +  EIxc  r’L^-f(hwf  +  GJ f  - jf'{hw f  -  crTenn2 

Cl  TC 


v_ 

77 1 


J717T  r  (  \  Cl  p"  /  \ 

- /(* 3 - /  (*3) 

a  ni7t 


~<2 


Tenn,  =  -?—  +  t— - 


+  2(\-v)[f(x3)f'(x3)  +  f(x3) 
b  r  „/  \2  ,  .  b~  »/_  \2  _  .  / —  ,  \  b~ 


>dx. 


2n  "  n 


hw  )2  +  2  A  [x  -  hf ) f  ( hw )  /'  ( hw )  + 1  f  f  ( hw ) 

71  J  71  77  77 


(56) 


where  the  omitted  limits  of  integration  are  from  x3  =  0  to  x3  =  hw . 
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Applying  the  stability  criteria  (31)  yields 


<7  <  ■ 


1  _  ,  ( mb  a  \ 
-Dp  —  +  — 

2  p  ^  a  mb , 


+  Term,  +  EI„  f  (h,  f  +  GJ ,  A f(hj 


tb 


2 n~  "  n 

b2 

TC 


P~  +  K  -T  J  /  (*3  fdx 3  +  Af  —  /  ( hw  f  +2Af(x-hf)-jf(  hw )  f'(hw  )  +  K  )2 


Tenn,  =  f 

77"  J 


/777T  r  (  \  O'  rH  /  \ 

— /  ( A3 ) - /  (*3) 

a  m/r 


“|2 


+  2(l-v)[f(x3)f'(x3)  +  f(x3y 
(57) 


\dx. 


where  a  proper  selection  of  /  (x3 )  minimizes  the  upper  bound  for  the  buckling  load. 

Notice  that  in  the  case  of  a  plate  panel  with  no  stiffener,  only  the  first  tenns  of  the 
numerator  and  denominator  in  (57)  remain  and  the  buckling  load  value  is 


<7  < 


'  Db 

2  p 


mb  a 

—  +  — 

v  a  mb 


J  _ 


2n2 


=  D. 


n 


2  f 


mb  a 

- + - 

v  a  mb 


y 

J 


(58) 


where  m  is  one  of  the  two  integers  closest  to  alb.  Expression  (58)  agrees  with  the  well- 
known  solution  for  the  buckling  of  a  plate  that  has  length  a ,  width  b ,  and  is  simply 
supported  at  all  edges. 

When  there  is  a  thin-webbed  stiffener  but  no  flange,  the  web  deflection  assumes 
the  form 


fM 


7TX, 

~b 


(59) 
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Figure  12  Plate  and  web  deflection  when  there  is  no  flange 

and  the  buckling  load  is  the  value 


3.  General  Buckling  Mode  and  Load 

Common  conjectures  for  /  (x, )  include 


i'  (  \  v  •  Tt'X's 

f  (x3)  =  — sin — - 
u  b  h 

w 

(61) 

r  (  \  ■ 

/2(x3)  =  sin  b 

(62) 

(63) 

Expressions  (61)  -  (63)  satisfy  boundary  conditions  (54)  and  (55). 

/,(0)=/2(0)=/J(0)  =  0 
k(0)=//(0)  =  //(0)  =  ^ 

When  (61)  -  (63)  are  substituted  into  (53),  they  yield  the  following  stiffener  deflection 
approximations: 
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Vl{XVXi) 

r2(xi,x3) 

k3(x,,x3) 


=  csin 


=  csin 


=  csin 


mnxx 

a 

mnxx 

a 

mnx j 
a 


hw  .  7tx 3 
—sin — 1 
b  h,„ 


sin- 


7TX, 


nx , 


X 

y 


Linear  combinations  of  (61)  -  (63)  are  often  used  to  achieve  better  accuracy  in 
approximating  the  stiffener  deflection.  The  deflection  mode  function  can  be  assumed  in 
the  form  of  a  series  of  functions  with  undetennined  coefficients.  If  the  series  of  functions 
is  complete,  by  increasing  the  number  of  terms  in  the  series  the  approximation  for  the 
buckling  load  converges  from  above  to  the  exact  buckling  load  value.  Functions  (61)  - 
(63)  belong  to  a  family  of  functions  whose  members  satisfy  the  boundary  conditions. 
Functions  belonging  to  such  a  family  are  of  the  form 


Note  that 


f{x3) 


M, 


sin- 


nx^ 

MK 


0 


(64) 


lim/(x3)=  lim 


M 


Mr  •  nx3 

— -sin- 


M 


nx. 


‘w  J 


Continuing  with  the  Rayleigh-Ritz  method,  substitution  of  the  assumed  web 
deflection  (64)  into  the  general  buckling  load  expression  (57)  yields 
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where  the  value  of  ,u  is  chosen  to  minimize  the  expression.  One  method  of  detennining 
ji  is  by  setting  equal  to  zero  the  ordinary  derivative  of  the  quotient  expression  (65)  with 
respect  to  //  and  solving  for  fj. .  The  value  can  also  be  obtained  graphically.  Three 
special  values  ju  =  1 ,  fi  -  2 ,  and  ju  — >  °°  relate  to  three  special  deflection  modes. 


For  /j  =  1 ,  the  deflection  mode  becomes  /  (x3 ) 


hw  .  nx^ 
—sin — 1 
b  K, 


and  the  buckling 


load  is 


cr  < 
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f  a 
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(66) 


Figure  13  depicts  the  deflection  of  the  structure  corresponding  to  mode  1  that  is  best 
described  by  a  bending  of  the  web  in  one  half-wave  and  a  rotation  of  the  flange  about  a 
point  at  the  top  of  the  web.  There  is  no  flange  bending  along  its  length.  Such  a  mode 
could  occur  when  the  flange  bending  stiffness  is  large  compared  to  flange  torsion  and 
web  bending.  An  example  of  this  occurrence  is  a  wide  T  flange  combined  with  a  thin 
web  possessing  relatively  low  flexural  rigidity. 


Figure  13  Mode  1  deflection 
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For  n  =  2 ,  the  deflection  mode  becomes  /  (x, ) 
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Figure  14  depicts  the  deflection  of  the  structure  corresponding  to  mode  2.  This  deflection 
corresponds  to  a  significant  bending  of  the  flange  and  web  with  no  flange  torsion.  This 
mode  could  occur  when  the  flange  torsional  stiffness  is  large  compared  to  flange  and  web 
bending. 


Figure  14  Mode  2  deflection 
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Figure  15  Mode  3  deflection 

Figure  15  depicts  the  deflection  of  the  structure  corresponding  to  mode  3.  This  deflection 
corresponds  to  the  flange  exhibiting  a  combination  of  bending  and  twisting  while  the  web 
tends  to  remain  straight.  Such  a  mode  could  occur  when  the  web  flexural  stiffness  is 
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large  compared  to  flange  bending  and  torsion  effects.  This  case  is  likely  to  occur  when 
the  flange  offers  little  or  no  additional  stiffness  to  the  plate  structure  compared  to  the 
contribution  of  the  web. 

The  value  of  /u  can  be  interpreted  to  indicate  the  deflection  behavior  based  on  the 
three  special  deflection  modes. 

Mode  ju  Flange  bending  Flange  torsion  Web  bending 


1 

1 

No 

Yes 

Yes 

2 

2 

Yes 

No 

Yes 

3 

oo 

Yes 

Yes 

No 
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4.  Buckling  Expression  Summary 

The  following  summarizes  the  analytic  expressions  that  predict  the  critical 
buckling  load  and  buckling  behavior  of  a  rectangular  stiffened  plate  with  a  thin  webbed 
stiffener  and  solid  flange  treated  as  a  beam.  The  general  expression  involves  a  constant 
//  determined  graphically.  The  value  of  ,u  indicates  the  deflection  behavior. 


General  Buckling  Load  Expression  (//  ^  0) 
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Mode  3  Expression  (//->“) 
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The  three  special  mode  expressions  provide  initial  approximations  for  the  upper  bound  of 

the  buckling  load  value.  Once  the  proper  value  of  //  is  estimated,  the  general  expression 

provides  a  value  that  is  less  than  or  equal  to  the  lowest  initial  approximation. 
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5.  Example 

As  an  example,  consider  a  stiffened  plate  structure,  simply  supported  on  all  sides, 
subjected  to  an  axial  compression  load  along  the  longitudinal  axis,  with  the  following 
parameters. 

Material  properties: 

E  =  30, 000  ksi  v  =  0.3 

Plate  properties: 

a  =  12  inches  =1829  mm 
b  =  20  inches  =  508  mm 
tp  =  0.3125  inches  =  7.9375  mm 

Web  properties: 

hs  =  4.1338  inches  =  105  mm 
tw  =  0.2362  inches  =  6  mm 

Bulb-flat  flange  properties: 

tbf  =  0.6693  inches  =  17  mm 
r  =  0.1969  inches  =  5  mm 
rx  =0.0236  inches  =0.6  mm 
a=n!6  radians  =  30° 


A  Maple  8  worksheet  provides: 


Dp  =  83840 
hf  =0.8106 
4  =  0-0309 
EIxc  =  926071 


Dw  =36212 
Af  =  0.5202 
J,  =0.0321 
GJf  =  370796 


m  =  4 
x  =  0.3214 
If  =0.2100 


K  =  3.3233 
y  =  0.3694 
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The  following  figure  shows  the  graph  of  the  general  expression  for  the  buckling 
load  of  the  plate  with  a  bulb-flat  stiffener. 


Figure  16  Maple  8  graph  of  the  buckling  load  for  the  bulb-flat 
stiffened  plate  example 
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Table  6  Results  for  the  bulb-flat  stiffened  plate  example 


Description 

a  <  Buckling  Load 

(psi) 

Remarks 

Plate  without  stiffener 

26773 

(58) 

Plate  stiffened  with  no  flange 

27581 

(60) 

3.02%  above  the  no  stiffener 
value. 

Mode  1 

44380 

(66) 

Mode  2 

31262 

(67) 

Mode  3 

31801 

(68) 

General  Expression 

27802 

Graphically  and  (65). 

jU  =  3.02 

3.84%  above  the  no  stiffener 
value. 

0.80%  above  the  no  flange  value. 

Without  stiffening,  the  plate  is  predicted  to  buckle  at  a  load  less  than  or  equal  to  26773 
psi.  When  a  flangeless  web  of  height  hw  =3.3233  in.  is  added  to  the  plate  panel,  the 
plate  buckles  at  a  load  of  less  than  or  equal  to  27581  psi.  Based  on  the  value  for  [i ,  the 
buckling  behavior  is  predicted  to  exhibit  flange  bending  and  torsion  with  relatively  little 
web  bending.  The  given  flange  cross-section  adds  relatively  little  additional  stiffness 
compared  to  the  web  stiffness. 

E.  FINITE  ELEMENT  ANALYSIS 

Numerical  and  experimental  investigations  give  insight  about  the  behavior  of 
stiffened  plates  and  serve  as  a  basis  of  comparison  to  theoretical  predictions.  The  results 
of  finite  element  analyses  for  several  stiffened  plate  models  using  MSC  Nastran 
2001/Patran  2001-r3  software  are  compared  to  analytical  predictions  to  establish  a  degree 
of  confirmation. 
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Figure  17  Stiffened  plate  pre-buckling  and  buckling  mode 

Several  simply  supported  rectangular  plate  structures  serve  as  finite  element 
models.  Each  plate  structure  has  a  plate  panel  of  length  a  ,  width  2b  ,  and  thickness  t 

such  that  t  <a  .  The  stiffener  divides  the  width  of  the  plate  in  halves.  Each  plate 
panel  half  extends  b  units  of  length  from  the  web  to  the  simply  supported  edge.  The 
stiffener  includes  a  thin  web  of  length  a ,  height  hw ,  and  thickness  tw  such  that 
tw  <sc  hw  <  a  .  Four  flange  configurations  are  examined:  no  flange,  bulb-flat  flange, 

circular  flange,  and  T-flange.  Each  stiffened  plate  structure  is  subjected  to  axial 
compression  due  to  a  uniform  normal  stress  <7,  which  increases  to  the  buckling  load. 
The  boundary  of  the  plate  panel  is  simply  supported  at  all  edges.  The  stiffener  is  simply 
supported  at  the  ends. 

The  derived  analytical  expressions  in  the  previous  section  were  modeled 
assuming  a  plate  panel  of  width  b  .  This  coincides  with  the  analysis  of  a  single  plate  unit 
containing  a  single  longitudinal  stiffener  that  would  be  part  of  a  larger  plate  structure 
where  b  units  of  length  separate  each  stiffener.  In  this  discussion,  the  plate  structure  is 
modeled  with  width  2b  ,  as  shown  in  Figure  17.  The  boundary  conditions  are  now 

W  (xj  ,-b)  =  W  (xj  ,b)  =  0 

W22  (XV-h)  =  W22(xi,b)  =  0 
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Consistency  with  the  finite  element  model  analysis  of  the  2b  width  model  requires 
changing  the  limits  of  integration  associated  with  the  plate  energy  terms  C'plalc  and  7/late 

from  x2=±b  /  2  to  x2=±b .  The  fundamental  problem  and  its  mode  shapes  remain 
unchanged. 

A  secondary  aim  of  the  finite  element  analysis  is  to  estimate  the  effects  of  several 
variables  or  factors.  The  intention  is  to  discern  the  sensitivity  of  the  analytic  expression 
accuracy  due  to  the  variability  of  web  height,  flange  area,  and  flange  configuration.  The 
web  height  variable  is  examined  at  two  levels,  the  flange  area  at  four  levels,  and  the 
flange  configuration  at  four  levels.  The  plate  panel  parameters  for  all  models  remain 
fixed  at  the  following  values: 

a  =  72  inches  =1829  mm 
b  =  20  inches  =  508  mm 
tp  =  0.3125  inches  =  7.9375  mm 

A  105  mm  stiffener  with  a  bulb-flat  flange  (see  model  1  from  Table  2)  has  a  web  height 
of  84.41  mm.  This  is  the  web  height  used  for  the  first  set  of  plate  models.  A  120  mm 
stiffener  with  a  bulb-flat  flange  (see  model  1  from  Table  2)  has  a  web  height  of  99.41 
mm.  This  is  the  web  height  used  for  the  second  set  of  plate  models.  The  web  thickness 
is  fixed  at  6  mm. 

Web  height  1:  hw  =  84.41  mm 
Web  height  2:  hw  =  99.41  mm 
Web  thickness:  tw  =  6  mm 

There  are  four  flange  areas  examined.  The  first  case  is  that  for  a  stiffener  with  no  flange. 
The  remaining  flange  area  values  are  derived  from  the  following  bulb-flat  flange 
parameters  (units  in  mm)  where  r\  =  (1/  10)tw  and  a -n / 6 . 


Flange  Area  Level 

K 

tbf 

r 

Area  (mm2 ) 

1 

N/A 

N/A 

N/A 

0 

2 

6 

17 

5 

335.59 

3 

6 

20 

5 

401.69 

4 

6 

25 

5 

523.39 
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The  flange  parameters  used  for  area  levels  3  and  4  are  not  from  the  data  in  Table  2  or 
Table  3.  Analyses  of  these  bulb-flat  flanges  provide  an  assessment  of  the  prediction 
accuracy  for  non-standard  parameters.  The  four  flange  configurations  examined  are 

Flange  Configuration  Label  Description 

1  N/A  No  flange 

2  b  Bulb-flat 

3  c  Circular 

4  t  T-flange 

For  the  stiffened  plate  structure  without  a  flange,  the  web  deflection  mode  is  of  the  form 
(59)  and  the  predicted  buckling  load  is  detennined  by 


DPb 


a . < 


( mb  a 
- + - 

v  a  mb 


+  Dh 


mnh 

vt 

y  a  ) 


+  2(1-1') 


n 2  3 


Otherwise  the  stiffener  has  a  flange  and  the  web  deflection  is  of  the  form  (64)  with  the 
predicted  buckling  load  detennined  by 
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Using  the  values  from  the  previous  two  pages  and  converting  the  units  from 
millimeter  to  inches  provides  the  parameters  for  the  stiffened  plate  models.  The 
following  tables  present  the  converted  parameter  values.  The  data  tables  in  Appendix  C 
for  the  stiffened  plate  finite  element  models  are  in  inch  units  (See  Table  40  -  Table  55). 
Converting  the  results  back  into  millimeter  units  may  yield  values  that  differ  slightly 
from  those  presented  in  previous  sections. 


Table  7  Summary  of  stiffened  plate  model  1  parameters 


Configuration 

Flange  Parameters  and  Properties 

Web  Parameters 

Model  1 

Af  =  0.5202  in2 

hw  =  3.3233  in 
tw  =  0.2362  in 

No  flange 

See  Table  40 

No  flange  area 

Same  as  above 

Bulb-flat  (b) 

See  Table  427 

tw  =  0.2362  in 
thf  =  0.6693  in 

r  =  0.1969  in 

i\  =(1/1  OX,.  =  0.0236  in 

a  =  7t  1 6  radians 
hf  =  0.8106  in 

hs  =4.1339  in 

Circular  (c) 

See  Table  48 

radius  =  0.4069  in 
hf  =  0.8138  in 

hs  =4.1371  in 

T-flange  (t) 

See  Table  54 

width  =  2.2020  in 

thickness  =  0.2362  in 

hs  =  3.5595  in 

7  The  parameters  for  this  model  configuration  were  used  in  a  grillage  tested  by  the  Naval  Surface 
Warfare  Center.  The  flange  and  web  parameters  are  the  same  as  used  in  the  example  on  page  52. 


58 


Table  8  Summary  of  stiffened  plate  model  2  parameters 


Configuration 

Flange  Parameters 

Web  Parameters 

Model  2 

Af  =  0.6226  in2 

hw  =  3.3233  in 
tw  =  0.2362  in 

No  flange 

See  Table  40 

No  flange  area 

Same  as  above 

Bulb-flat  (b) 

See  Table  43 

tw  =  0.2362  in 
thf  =  0.7874  in 

r  -  0.1969  in 

i\  =(1/10)6,.  =  0.0236  in 

a  =  n  1 6  radians 
hf  =  0.8788  in 

hs  =4.2021  in 

Circular  (c) 

See  Table  49 

radius  =  0.4452  in 
hf  =  0.8904  in 

hs  =4.2137  in 

T-flange  (t) 

See  Table  54 

width  =  2.6357  in 

thickness  =  0.2362  in 

h s  =  3.5595  in 

Table  9  Summary  of  stiffened  plate  model  3  parameters 


Configuration 

Flange  Parameters 

Web  Parameters 

Model  3 

Af  =0.8113  in2 

hw  =  3.3233  in 
tw  =  0.2362  in 

No  flange 

See  Table  40 

No  flange  area 

Same  as  above 

Bulb-flat  (b) 

See  Table  44 

tw  =  0.2362  in 
thf  =  0.9843  in 

r  =  0.1969  in 

rx  =(1/1  OX,.  =  0.0236  in 

a  =  7t  /  6  radians 
hf  =  0.9924  in 

hs  =4.3157  in 

Circular  (c) 

See  Table  50 

radius  =  0.5082  in 
hf  =  1.0163  in 

hs  =  4.3396  in 

T-flange  (t) 

See  Table  54 

width  =  3.4343  in  =  87.23  mm 

thickness  =  0.2362  in  =  6  mm 

h ?  =  3.5595  in 
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Table  10  Summary  of  stiffened  plate  model  4  parameters 


Configuration 

Flange  Parameters 

Web  Parameters 

Model  4 

Same  as  model  1 

hw  =  3.9138  in 
tw  =  0.2362  in 

No  flange 

See  Table  41 

No  flange  area 

Same  as  above 

Bulb-flat  (b) 

See  Table  45 

Same  as  model  lb 

h s  =  4.7244  in 

Circular  (c) 

See  Table  51 

Same  as  model  lc 

hs  =  4.7276  in 

T-flange  (t) 

See  Table  55 

Same  as  model  It 

h s  =4.1501  in 

Table  11  Summary  of  stiffened  plate  model  5  parameters 


Configuration 

Flange  Parameters 

Web  Parameters 

Model  5 

Same  as  model  2 

hw  =3.9138  in 
tw  =  0.2362  in 

No  flange 

See  Table  41 

No  flange  area 

Same  as  above 

Bulb -flat  (b) 

See  Table  46 

Same  as  model  2b 

hs  =  4.7926  in 

Circular  (c) 

See  Table  52 

Same  as  model  2c 

hs  =  4.8042  in 

T-flange  (t) 

See  Table  55 

Same  as  model  2t 

hs  =  4.1501  in 

Table  12  Summary  of  stiffened  plate  model  6  parameters 


Configuration 

Flange  Parameters 

Web  Parameters 

Model  6 

Same  as  model  3 

hw  =3.9138  in 
tw  =  0.2362  in 

No  flange 

See  Table  41 

No  flange  area 

Same  as  above 

Bulb-flat  (b) 

See  Table  47 

Same  as  model  3b 

hs  =  4.9063  in 

Circular  (c) 

See  Table  53 

Same  as  model  3c 

hs  =4.9301  in 

T-flange  (t) 

See  Table  55 

Same  as  model  3t 

hs  =  4.1501  in 
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The  following  tables  present  intermediate  results  for  the  stiffened  plate  models. 


Table  13  Model  1  intermediate  results 


Result 

Bulb-flat 

Circular 

T-flange 

hf  (in) 

0.8106 

0.8138 

0.2362 

x  (in) 

0.3214 

0.4069 

0.1181 

y  (in) 

0.3694 

0.4069 

1.1010 

/f(in4) 

0.0526 

0.0431 

0.2126 

4  (in4) 

0.0309 

0.0215 

0.2102 

J f  (in4) 

0.0321 

0.0431 

0.0090 

If  (in4) 

0.2100 

0.1292 

0.2199 

EI*c  (lb-in2) 

926071 

645945 

6305666 

GJ f  (lb-in2) 

370796 

496881 

104091 

Mode  1  (psi) 

35599 

36107 

34548 

Mode  2  (psi) 

29053 

28905 

31902 

Mode  3  (psi) 

29356 

29455 

35105 

Table  14  Model  2  intermediate  results 


Result 

Bulb-flat 

Circular 

T-flange 

hf  (in) 

0.8788 

0.8904 

0.2362 

x  (in) 

0.3413 

0.4452 

0.1181 

y  (in) 

0.4145 

0.4452 

1.3179 

I A  in4) 

0.0773 

0.0617 

0.3633 

4  (in4) 

0.0473 

0.0308 

0.3604 

J f  (in4) 

0.0449 

0.0617 

0.0109 

If  (in4) 

0.3119 

0.1851 

0.3720 

EIXC  (lb-in2) 

1417646 

925443 

10813397 

GJ  A  lb-in2) 

517932 

711879 

126078 

Mode  1  (psi) 

36164 

36946 

34614 

Mode  2  (psi) 

29261 

29001 

34229 

Mode  3  (psi) 

30474 

30553 

40794 
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Table  15  Model  3  intermediate  results 


Result 

Bulb-flat 

Circular 

T-flange 

hf  (in) 

0.9924 

1.0163 

0.2362 

x  (in) 

0.3753 

0.5082 

0.1181 

y  (in) 

0.4886 

0.5082 

1.7172 

/f(in4) 

0.1361 

0.1047 

0.8011 

4  (m4) 

0.0872 

0.0524 

0.7974 

J f  (in4) 

0.0735 

0.1047 

0.0144 

If  (in4) 

0.5564 

0.3142 

0.8125 

EIxc  (lb-in2) 

2617311 

1571197 

23921250 

a./,  (lb-in2) 

848126 

1208613 

166564 

Mode  1  (psi) 

37428 

38882 

34713 

Mode  2  (psi) 

29798 

29247 

41024 

Mode  3  (psi) 

33128 

33117 

57198 

Table  16  Model  4  intermediate  results 


Result 

Bulb-flat 

Circular 

T-flange 

hf  (in) 

0.8106 

0.8138 

0.2362 

x  (in) 

0.3214 

0.4069 

0.1181 

y  (in) 

0.3694 

0.4069 

1.1010 

/f(in4) 

0.0526 

0.0431 

0.2126 

4  (in4) 

0.0309 

0.0215 

0.2102 

Jf  (in4) 

0.0321 

0.0431 

0.0090 

If  (in4) 

0.2100 

0.1292 

0.2199 

^4  (lb-in2) 

926071 

645945 

6305666 

GJf  (lb-in2) 

370796 

496881 

104091 

Mode  1  (psi) 

34602 

35109 

33551 

Mode  2  (psi) 

28862 

28657 

32783 

Mode  3  (psi) 

29541 

29495 

37813 
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Table  17  Model  5  intermediate  results 


Result 

Bulb-flat 

Circular 

T-flange 

hf  (in) 

0.8788 

0.8904 

0.2362 

x  (in) 

0.3413 

0.4452 

0.1181 

y  (in) 

0.4145 

0.4452 

1.3179 

/c(in4) 

0.0773 

0.0617 

0.3633 

4  (in4) 

0.0473 

0.0308 

0.3604 

J f  (in4) 

0.0449 

0.0617 

0.0109 

I  f  (in4) 

0.3119 

0.1851 

0.3720 

EIXC  (lb-in2) 

1417646 

925443 

10813397 

GJ  f  (lb-in2) 

517932 

711879 

126078 

Mode  1  (psi) 

35166 

35947 

33618 

Mode  2  (psi) 

29148 

28791 

35980 

Mode  3  (psi) 

30844 

30672 

45530 

Table  18  Model  6  intermediate  results 


Result 

Bulb-flat 

Circular 

T-flange 

hf  (in) 

0.9924 

1.0163 

0.2362 

x  (in) 

0.3753 

0.5082 

0.1181 

y  (in) 

0.4886 

0.5082 

1.7172 

/c(in4) 

0.1361 

0.1047 

0.8011 

4  (in4) 

0.0872 

0.0524 

0.7974 

J f  (in4) 

0.0735 

0.1047 

0.0144 

If  (in4) 

0.5564 

0.3142 

0.8125 

^4  (lb-in2) 

2617311 

1571197 

23921250 

GJ  f  (lb-in2) 

848126 

1208613 

166564 

Mode  1  (psi) 

36430 

37882 

33718 

Mode  2  (psi) 

29886 

29129 

45307 

Mode  3  (psi) 

33964 

33430 

67754 
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The  following  pictures  show  the  deflection  of  the  plate  and  web  at  the  onset  of 
buckling  for  Model  lb.  The  plate  and  web  panels  buckle  in  4  half  waves  along  the 
longitudinal  axis.  The  bulb-flat  flange  is  not  visible,  but  its  properties  exist  in  the  curve 
that  defines  the  top  of  the  web. 


Figure  18  MSC  Nastran  picture  1  for  stiffened  plate  Model  lb 


Figure  19  MSC  Nastran  picture  2  for  stiffened  plate  Model  lb 
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Circular 

T-flange 


Figure  20  Maple  8  graph  1  of  the  buckling  loads  for  Model  1 
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Figure  21  Maple  8  graph  2  of  the  buckling  loads  for  Model  1 
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The  pictures  below  show  the  difference  in  the  buckling  behavior  between  models 
2b  and  3b.  In  each  picture,  the  top  of  the  web  has  the  beam  characteristics  of  a  bulb-flat 
flange.  The  flange  in  the  left  picture  has  an  area  property  of  0.6226  in  .  The  flange  in  the 


Figure  23  MSC  Nastran  pictures  showing  the  buckling  behavior  of  the  webs 
with  bulb-flat  flanges  for  Models  2b  (left)  and  3b  (right) 

right  picture  has  an  area  property  of  0.8113  in  .  These  properties  are  not  visible  in  the 
pictures,  but  exist  in  the  curve  that  defines  the  top  of  each  web.  The  web  height  in  both 
pictures  is  3.3233  in.  From  Table  20,  the  //  values  for  the  left  and  right  pictures  are  2.57 
and  2.25  respectively.  This  indicates  there  is  more  lateral  bending  of  the  Model  2b  flange 
at  the  onset  of  buckling,  which  agrees  with  the  pictures  shown  above. 
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The  pictures  below  show  the  difference  in  the  buckling  behavior  between  models 

4b,  5b,  and  6b  from  left  to  right  respectively.  The  flanges  have  area  properties  of  0.5202 

2  2  2 

in  ,  0.6226  in“,  and  0.81 13  in  ,  respectively.  The  web  height  in  the  pictures  is  3.9138  in. 
From  Table  20,  the  //  values  are  2.68,  2.40,  and  2.11,  respectively. 


Figure  24  MSC  Nastran  pictures  showing  the  buckling  behavior  of 
the  webs  with  bulb-flat  flanges  for  Models  4b  (left),  5b  (center),  and  6b 
(right) 


Figure  25  MSC  Nastran  pictures  showing 
the  buckling  behavior  of  the  webs  with  T-flanges 
for  Models  4t  (left)  and  6t  (right) 


The  pictures  above  show  the  difference  in  the  buckling  behavior  between  models  4t  and 
6t.  The  flanges  have  area  properties  of  0.5202  in"  and  0.81 13  in  ,  respectively.  The  web 
height  in  the  pictures  is  3.9138  in.  From  Table  22,  the  //  values  are  1.51  and  1.08, 
respectively. 
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The  following  tables  show  the  formula  predictions  and  finite  element  results. 


Table  19  No  flange  FE  results 


Model# 

FE  Results  (psi) 

Analytic  Prediction  (psi) 

Percent  Change 

1-3 

26928 

27181 

0.94% 

4-6 

26892 

27138 

0.91% 

Table  20  Bulb-flat  FE  results 


Model# 

FE  Results  (psi) 

Analytic  Prediction  (psi) 

A 

Percent  Change 

lb 

27672 

28397 

2.86 

2.62% 

2b 

27877 

28824 

2.57 

3.40% 

3b 

28203 

29642 

2.25 

5.10% 

4b 

27663 

28425 

2.68 

2.75% 

5b 

27881 

28908 

2.40 

3.68% 

6b 

28217 

29852 

2.11 

5.79% 

Table  21  Circular  FE  results 


Model# 

FE  Results  (psi) 

Analytic  Prediction  (psi) 

A 

Percent  Change 

lc 

27558 

28246 

2.78 

2.50% 

2c 

27703 

28524 

2.54 

2.96% 

3c 

27936 

29000 

2.29 

3.81% 

4c 

27498 

28178 

2.65 

2.47% 

5c 

27654 

28469 

2.44 

2.95% 

6c 

27907 

28996 

2.21 

3.90% 

Table  22  T-flange  FE  results 


Model# 

FE  Results  (psi) 

Analytic  Prediction  (psi) 

A 

Percent  Change 

It 

28767 

31864 

1.85 

10.77% 

2t 

29030 

33200 

1.45 

14.36% 

3t 

29291 

34220 

1.15 

16.83% 

4t 

28739 

32146 

1.51 

11.85% 

5t 

28933 

32945 

1.23 

13.87% 

6t 

29119 

33460 

1.08 

14.91% 
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Table  19  compares  the  results  for  the  stiffened  plates  without  flanges.  The 
formula  values  nearly  agree  with  the  finite  element  results.  Table  20  compares  the  bulb- 
flat  configurations  for  each  model.  The  formula  values  are  less  than  6%  above  the  finite 
element  results.  Examination  of  the  ,u  values  indicates  that  the  buckling  behavior  of  the 
stiffened  plates  tend  to  have  significant  bending  of  the  flange  and  web  with  a  degree  of 
flange  torsion.  This  type  of  buckling  behavior  can  be  seen  in  the  Nastran  pictures  on  the 
previous  pages.  As  the  flange  areas  increase  between  models  1,  2,  and  3,  and  also 
between  models  4,  5,  and  6,  the  //  values  tend  to  approach  the  value  2,  indicating 
perhaps  an  increasing  torsional  rigidity  effect.  The  buckling  loads  tend  to  increase  as  the 
flange  area  increases,  as  expected.  Overall,  the  results  indicate  that  at  the  onset  of 
buckling,  the  bulb-flat  stiffened  plate  structures  deflect  with  a  combination  of  flange 
bending  and  torsion  along  with  web  bending.  Table  21  compares  the  circular  flange 
configurations,  which  appear  to  buckle  at  a  slightly  lower  value  than  the  bulb-flat 
configurations,  but  behave  similar  to  them  at  the  onset  of  buckling.  Table  22  compares 
the  T-flange  configurations.  From  the  finite  element  results,  the  buckling  load  of  a  plate 
with  a  bulb-flat  stiffener  is  3%  -  4%  less  than  that  of  a  plate  with  a  T-flange  stiffener  with 
the  same  cross-sectional  area.  Based  on  the  fu  values  for  the  T-flange  configuration, 
bending  of  the  flange  tends  to  decrease  as  the  cross-sectional  area  increases.  The 
difference  in  the  formula  values  from  the  finite  element  values  could  be  due  to 
fundamental  assumption  (ix).  Though  the  T-flange  is  solid,  treating  the  flange  as  a  thin 
web  plate  strip  instead  of  a  beam  may  achieve  more  accurate  predictions.  This  may 
explain  the  increase  in  error  for  the  bulb-flat  and  T-flange  cases  as  the  area  increases  due 
to  flange  widening.  The  increasing  error  is  less  pronounced  for  the  circular  flange 
models. 
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IV.  FINDINGS  AND  CONCLUSIONS 


The  following  summarizes  the  major  findings  and  conclusions  of  this  work.  Each 
major  finding  or  conclusion  is  stated  followed  by  explanatory  comments. 

•  Determined  cross-sectional  boundary  equations  for  the  bulb-flat  flange  cross- 
section. 

In  order  to  conduct  an  accurate  analysis  of  the  bulb-flat  cross-section,  it  is 
necessary  to  determine  equations  that  define  or  at  least  approximate  the  boundary  that  is  a 
closed  bounded  plane  region.  The  equation  that  defines  the  boundary  of  a  simple 
geometric  cross-section  like  a  circle  or  ellipse  is  well  known.  The  boundary  equation  for 
uncommon  and  asymmetric  cross-sections  is  often  very  difficult  to  define.  As  a  result  of 
this  study,  functions  exist  that  define  the  boundary  of  the  bulb-flat  flange  cross-section. 

•  Derived  planar  property  value  expressions. 

Determining  the  critical  buckling  load  of  a  stiffened  plate  structure  requires 
knowledge  of  several  planar  property  values  of  the  plate’s  cross-section.  Based  on  the 
boundary  equations  and  the  application  of  multivariable  calculus,  the  double  integral 
provides  expressions  for  calculating  the  various  planar  properties  of  the  bulb-flat  flange 
cross-section.  Though  the  integral  expressions  are  complicated,  they  could  be  simplified 
by  assigning  fixed  values  to  certain  variables. 

•  Determined  an  approximate  torsional  constant  expression  that  is  more 
accurate  than  idealizing. 

The  Saint- Venant  torsional  constant  is  one  of  the  key  property  values  involved  in 
the  analysis  of  stiffened  plate  structures.  Due  to  the  uncommon  shape  and  asymmetrical 
property  of  the  bulb-flat  cross-section,  determining  the  exact  expression  for  the  Saint- 
Venant  torsional  constant  is  difficult  and  does  not  exist  in  published  literature.  As  a 
result  of  this  study,  an  approximate  expression  exists  for  the  Saint- Venant  torsional 
constant  of  a  specified  class  of  bulb-flat  cross-sections.  The  approximate  expression  is  a 
one -term  function  relating  the  cross-sectional  area  to  the  polar  moment  of  inertia.  The 
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approximate  expression  is  more  accurate  than  the  estimate  obtained  by  idealizing  the 
cross-section  as  an  angle  flange. 

•  Demonstrated  that  the  torsional  property  of  the  bulb-flat  stiffener  is  better  than 
previously  understood.  The  torque-carrying  capacity  of  a  bulb-flat  stiffener 
(possessing  no  structural  flaws)  is  greater  than  that  of  an  area-equivalent  angle 
stiffener. 

In  previous  investigations  the  torsional  nature  of  the  bulb-flat  stiffener  was 
idealized  as  an  angle  flange  stiffener.  That  is,  the  bulb-flat  flange  cross-section  was 
treated  like  a  rectangular  cross-section  in  regards  to  the  torsional  and  warping  properties. 
Such  treatment  imputed  error  in  the  calculation  of  the  bulb-flat  stiffener’s  torsional 
rigidity  resulting  in  conservative  estimates.  Finite  element  analysis  indicates  that  the 
torque-carrying  capacity  of  a  bulb-flat  stiffener  (possessing  no  structural  flaws)  is  greater 
than  that  of  an  area-equivalent  angle  stiffener. 

•  Derived  a  general  expression  to  predict  the  buckling  load  due  to  the  stiffener 
tripping  of  a  simply  supported  rectangular  stiffened  plate  subjected  to  axial 
compression.  In  the  investigation,  the  predicted  value  is  less  than  6%  higher 
than  the  finite  element  result. 

Understanding  the  elastic  stability  of  stiffened  plate  structures  is  important  to  the 
analyst,  designer,  and  educator  involved  in  the  analysis  and  design  of  structures.  Use  of 
the  energy  method  provides  a  technique  to  derive  a  general  expression  for  the  buckling 
load  due  to  the  stiffener  tripping  of  a  simply  supported  rectangular  stiffened  plate 
subjected  to  axial  compression.  As  a  result  of  this  study,  a  useful  analytic  expression 
exists  that  allows  the  user  to  predict  the  critical  buckling  load  and  the  buckling  behavior 
of  a  stiffened  plate  at  the  onset  of  stiffener  tripping.  The  onset  of  stiffener  tripping 
negates  the  stiffener’s  support  to  the  plate  panel  and  leads  to  eventual  collapse  of  the 
structure.  The  general  expression  involves  a  constant  called  ,u  that  is  determined 
graphically.  The  value  of  ,u  can  be  interpreted  to  indicate  the  deflection  behavior  based 
on  three  special  deflection  modes. 
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Mode 

Flange  bending 

Flange  torsion 

Web  bending 

1 

1 

No 

Yes 

Yes 

2 

2 

Yes 

No 

Yes 

3 

oo 

Yes 

Yes 

No 

Finite  element  analyses  of  several  stiffened  plate  models  is  presented  to  validate 
the  analytic  expressions  and  assess  the  sensitivity  of  the  formulas  to  factor  variability. 
Knowledge  of  MSC  Nastran/Patran  software  and  finite  element  principles  along  with  the 
Maple  8  computer  environment  was  essential  for  this  investigation. 

•  Demonstrated  that  the  buckling  behavior  of  the  bulb-flat  stiffened  plate  is 
markedly  unlike  that  of  the  T-flange  configuration.  The  bulb-flat  tends  to 
buckle  more  laterally  and  have  a  lower  buckling  load  value  than  an  area- 
equivalent  T-flange  stiffened  plate. 

The  bulb-flat  stiffened  plate  tends  to  buckle  in  a  Mode  2  fashion,  while  the  T- 
flange  stiffened  plate  tends  to  buckle  in  a  Mode  1  fashion.  Hence,  the  torsional  superior 
bulb-flat  tends  to  bend  laterally  at  the  onset  of  stiffener  tripping,  while  the  flexural 
superior  T-flange  tends  to  twist  axially  at  the  onset  of  stiffener  tripping. 

As  a  result  of  this  investigation,  several  essential  future  research  directions  exist. 
Future  research  should  include: 

•  Investigating  methods  to  determine  the  ,u  value  by  other  than  graphical 
means. 

•  Investigating  the  use  of  conformal  mapping  to  determine  the  exact  expression 
for  the  bulb-flat  torsional  constant. 

•  Investigating  solutions  to  the  torsion  problem  for  asymmetric  cross-sections. 

•  Conducting  investigations  of  other  flange  cross-sections,  multiple  stiffener 
configurations,  and  grillages. 

•  Conducting  further  investigations  of  the  T  flange  by  treating  the  flange  as  a 
thin  web  plate  strip  instead  of  a  beam.  This  treatment  may  achieve  more 
accurate  predictions  for  the  T-flange  buckling  loads. 

•  Developing  algorithms  that  increase  efficiency  in  stiffened  plate  analysis  and 
design. 
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APPENDIX  A.  THE  BULB-FLAT  FLANGE  HEIGHT  FORMULA 


This  appendix  presents  the  derivation  of  the  formula  for  the  height  of  the  bulb-flat 
flange  cross-section.  Maple  8  computer  environment  software  was  used  to  derive  some 
of  the  formula  final  forms.  Five  independent  variables  uniquely  detennine  a  specific 
bulb-flat  cross-section. 

tw  Thickness  of  the  web  (length) 

tbf  Thickness  of  the  flange  bulb  (length) 

r  Radius  of  curvature  of  flange  bulb  and  neck  curve  (length) 

rx  Radius  of  curvature  of  the  cross-section  comer  (length) 

a  Slope  angle  of  the  flange  neck  (radians) 

It  will  be  shown  that  only  three  of  the  independent  variables,  thf ,  r  and  a ,  determine  the 

height  of  the  flange.  Throughout  this  discussion,  it  is  assumed  ()<a  <n  1 2  .  The 
following  figures  show  a  model  of  the  bulb-flat  flange  oriented  horizontally  with  the  bulb 
extended  to  one  side  of  the  web. 
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Figure  26  Bulb-flat  flange  geometry 


Figure  27  MSC  Patran  Bulb-flat  flange  boundary 
points  for  cross-section  properties 
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The  three  circles  in  Figure  26  show  the  radii  of  curvature  at  the  corners  of  the  flange. 
The  height  of  the  flange,  denoted  h  f ,  is  the  distance  from  the  left  boundary  x  =  0  to  the 

right  boundary  x  =  h  f .  Consider  the  line  segment  tangent  to  the  two  circles  of  radius  r  . 

In  Figure  27,  the  line  segment  connects  point  7  (x7,y7 )  and  point  8  (x8,v8)  where 


and 


x7  is  to  be  determined 
y7  =  tw  +  r-r  sin  a 


xg  =  r  +  r  cos  a 
ys=tw  +  tbf-r  +  rsma 

The  flange  height  can  be  expressed  simply  as 

hf  =x7  +rcosa 

Once  x7  is  detennined,  hf  is  determined  from  (71).  The  slope  of  line  segment  is 


(69) 


(70) 


(71) 


,  Av  v7  -  Vo 
slope  =  —  =  — — —  =  -cota 


Ax  x7  -  xg 


*7  =-(t7  -y8)tan«  +  x8 


From  (69)  and  (70) 


From  (72)  and  (73) 


and 


T7  - Ts  =  (K  +  r- rsina) - [tw  +  thf  -r  +  r sina) 
=  -tbf  +  2r(l-sina) 


x7  -  x8  =  thf  tan  a  -  2 r  tan  a  +  2r  sec  a  -  2 r  cos  a 


x7  =  thf  tan  a  +  r-  2 r  tan  a  +  2 r  sec  a  -  r  cos  a 


Substituting  (75)  into  (71)  yields  the  following  fonnula. 

hf  =  hr  tana  +  r(l-  2  tan  a  +  2  sec  a) 


(72) 

(73) 

(74) 

(75) 

(76) 
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APPENDIX  B.  THE  PLANAR  PROPERTY  FORMULAS 


The  purpose  of  this  appendix  is  to  present  the  derivation  of  the  planar  property 
formulas  for  the  bulb-flat  flange.  It  is  assumed  0  <  a  <  TC !2  throughout  the  discussion. 

Let  the  closed  bounded  plane  region  R  be  defined  by  0  <  x  <  hf , 

g(x )  <  v  <  f(x),  with  / (x)  and  g(x)  on  x  e  [o,  hf  j  ,  such  that 


and 


where 


fix) 


fi  (x)  =  (fv  +  tbf~r)  +  Y  ~(x~rY  >  0  <x<r  +  r  cos  a 
f2  (x)  =  mx  +  b,  r  +  r  cos  a<x<hf-r  cos  a 


f3{x)  =  {tK  +  r)-^>'2  -(x~hf\  >  hf~r 


cos  a<x<hf 


g(x) 


gi{x)  =  '\  -{x-’\)2 »  °^x 

g2  (x)  =  0,  rx  <  x  <  hf 


<  r, 


m  =  -  cot  a 

b  =  tw  +  thf  +  r  (cot  a  +  esc  a- 1) 
hf  =  tbf  tana  +  r{\-  2  tan  a  +  2  sec  a) 


(77) 


(78) 


(79) 


By  the  theorem  of  integrability  of  a  piecewise  function8,  since  /  and  g  are  piecewise 
defined  functions  on  the  same  closed  interval,  then  /  and  g  are  integrable  on  the  same 
closed  interval.  The  planar  property  fonnulas  can  be  detennined  from  the  following 
expression  where  m,n  are  nonnegative  integers. 

jj x"y’flL4  =  j  J  xmyndydx  =  -^—  J  xm  |[/(x)]"+1 -[g(x)]"+1|  Jx 

R  0  g(*)  n  +  1  0 


8  Ross,  K.A.,  Elementary  Analysis:  The  Theory  of  Calculus,  Springer,  New  York,  pp.  258-259. 
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— -a 

\\xmyndA  =  —  J  (r  +  rsm9y  [tw  +tbf  -  r  +  rcosd^j  rcos0d0 

R  n  +  \  n 

2 

hf-rcosa 

If  /  \n+ 1 

+ -  xm  -x  cot  a  +  tw  +  tbf  +  /•  cot  a  +  r  esc  a  -  r )  dx 

n  + 1  ^ J  v 

r+rcosa 

i  o 

+ - -  J  [lif  -t-rsin#)  [tw  +  r-r  cos 0)"+' r  cos 0  d0 


- r\ 

n  + 1 


J  (l  +  sin^y"  (l-cos^),,+1  cos#  d0 


where 

hf  =  thf  tana  +  r(l-2tana  +  2seca) 

0  <a<~,  0  <x<hf 

2  2  2  f 

By  definition,  the  area  of  the  closed  bounded  plane  region  R  is  the  value  of  the  integral 

hr  /(*)  hr 

Area  =  jj  x°y°  dA  =  JJ  dA  =  J  J  dy  dx  =  J  / (x)-  g(x)  dx 

R  R  0  g(x)  0 

Hence, 

r+rcosa  h^-r  coscx  hj-  ^  h y 

JJ  dA  =  J  f  (x)dx+  J  f2  (x)dx+  J  /3  (x)dx-  J  gj  [x)dx-  J  g2  (x)dx 


Z  7 

JJfiL4=  J  (tw  +  tbf  -r  +  rcos#)rcos#  d0+  J  [~xcota  +  tw  +  tbf  +rcota  +  rcsca-r^dx 


+  (t.+r-r 


cos0)r  cos0  d0-r2  J  (l-cos#)cos 0  d0 


Integrating  and  collecting  terms  yields  the  area  of  the  bulb-flat  flange: 


Af=-  1 - [r2  +  r2 )  +  rtw  (l  -  2  taner  +  2secer)  +  rtbf  (l  -  taner  +  secct) 


1  2 

-I — tan  a  +  t,  t, ,  tan  a 
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Other  planar  property  formulas  are  determined  in  similar  fashion  and  are 
summarized  below  (algebraic  expansion  omitted). 


Mx  =  jj  xy'dA  =  JJ  y  dA 

R  R 

n 

— a 

1  2 

=  o  \  (o+Vr  +  r 


cos#)  r  cos  Odd 


+  —  J  ( -x  cot  a  +  tw  +  thf  +  r  cot  rz  +  r  esc  rz  -  r  j  dx 

^  r+rcosa 

1  °  1  ° 

+  —  J  [tw  +  r-rcos#)2  rcos#  dd —  r*  J  (l-cos#)2  cos#  dd 

^  71  ^  71 

a —  — 

2  2 

[y  -  JJ  xly°dA  =  JJ  x  dA 

R  R 

71 

— a 
2 

=  J  (r  +  rsmd)(tw+tbf  -r  +  rcosd^rcosd  dd 


|  x[-x  cot  a  +  tw  +  tbf  +  r  cot  a  +  r  esc  a  -  /•)  dx 


+  r  sin 0  ( t  +r  -r 


cos#)r  cos  6  d6-r*  J  (l  +  sin#)(l-cos#)cos#  dd 


83 


=  jjx°y2dA=jjy2  dA 

R  R 

71 

— a 

1  2p  /  \  3 

=  -  J  \tw  +  ty-r  +  rcos&)  rcos0d0 


+  —  J  (-xcota  +  tw  +  tbf  +  rcota  +  rcsca-rsj  dx 

^  r+rcosa 

1  °  1  ° 

+  t  J  {tw  +  r  -r cos#)3  r  cos#  dd  — 1\4  J  (l-cos#)3  cos#  d0 

^  n  ^  n 

a —  — 

2  2 

lyy  =  JJ  x2y°dA  =  JJ  X2  dA 

R  R 

71 

— a 
2 

=  J  (r  +  rsin0y  [tw  +  tbf  -r  +  rcos#)rcos#  d0 

71 

~2  ^ 
hr-rcosa 

+  J  x2  [~xcoia  +  tw  +tbf  +  rcota  +  rcsca-r'j  dx 

r+rcosa 

0  2  0 

+  J  (hf+rsin0)  (tw  +  r  -  r  cos  0)  r  cos  0  d0  -  rx  4  J  (l  +  sin#)"  (l  -cos0)cos0  d0 

71  71 

a —  — 

2  2 

Ky  =  JJ  x'y'dA  =  JJ XV  dA 

R  R 

71 

—~a 

=  —  J  {r  +  rsm0^(tw  +  thf  -r  +  rcos#)  rcos0d0 

2  jt 

d 

^  hj--rcosa 

+  —  J  x(-xcota  +  tw  +  tbf  +  rcota+rcsca-r^  dx  ( 

r+rcosa 

1  °  2 

+  -  J  (■ hf  +  rsin0)(tw  +  r  -r  cos #)"  rcos0  d0 

2  n 
a — 

2 

1  ° 

—  r [4  J  (l  +  sin#)(l  -cos#)  cos 0  d0 
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(89) 


Kyc  =  4,  -  Arxy  (90) 

Ic=Ixc+Iyc  (91) 
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APPENDIX  C.  DATA  TABLES 


This  appendix  presents  flange,  stiffener,  and  plate  data  tables  used  to  model  the 
cross-sections  for  finite  element  analysis.  When  creating  a  finite  element  model  of  a 
structure,  a  matrix  represents  the  structure.  The  information  in  the  left  portion  of  each 
table  list  parameters  or  references  related  to  the  set  of  points  contained  in  the  right 
portion  of  the  table. 

1.  DATA  FOR  BULB-FLAT  FLANGE  CROSS-SECTION 

See  Figure  3  and  Figure  5. 


Table  23  Data  for  Flange  Model  1  in  Table  2  (units  in  mm) 


K  =  6 

[0.0000  0.6000  0.0000][0.1757  0.1757  0.0000][0.6000  0.0000 

/  =  17 

0.0000] [20.5885  0.0000  0.0000][20.5885  6.0000  0.0000] [  1 8.0885 

bf  1  ' 

6.6699  0.0000][16.2583  8.5000  0.0000][9.3301  20.5000 

r  =  5 

0.0000][5.0000  23.0000  0.0000][0.0000  18.0000  0.0000] 

ii 

o 

s 

II 

On 

O 

a-n!6 

hf  =  20.59 

Af  =335.59 

Table  24  Data  for  Flange  Model  3  in  Table  2  (units  in  mm) 


t  =8 

w 

[0.0000  0.8000  0.0000][0.2343  0.2343  0.0000][0.8000  0.0000 

/  =  17 

0.0000][20.5885  0.0000  0.0000][20.5885  8.0000  0.0000] [  1 8.0885 

bf  1  ' 

8.6699  0.0000][16.2583  10.5000  0.0000] [9.330 1  22.5000 

r  =  5 

0.0000][5.0000  25.0000  0.0000][0.0000  20.0000  0.0000] 

o 

oc 

II 

* 

'K* 

S 

II 

5T 

a-nl6 

hf  =  20.59 

Af  =376.71 
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Table  25  Data  for  Flange  Model  5  in  Table  2  (units  in  mm) 


K  =  7 

[0.0000  0.7000  0.0000][0.2050  0.2050  0.0000][0.7000  0.0000 

t  =19 

0.0000][22.8205  0.0000  0.0000][22.8205  7.0000  0.0000][20.0705 

bf 

7.7369  0.0000][18.0574  9.7500  0.0000] [  1 0.263 1  23.2500 

r  =  5.5 

0.0000][5.5000  26.0000  0.0000][0.0000  20.5000  0.0000] 

II 

o 

*~~+ 

s 

II 

o 

a  =  nt  6 

hf  =  22.82 

Af  =422.19 

Table  26  Data  for  Flange  Model  7  in  Table  2  (units  in  mm) 


K  =  10 

[0.0000  1.0000  0.0000][0.2929  0.2929  0.0000][1.0000  0.0000 

t,r  =  19 

0.0000][22.8205  0.0000  0.0000][22.8205  10.0000 

bf  Ly 

0.0000][20.0705  10.7369  0.0000][18.0574  12.7500 

r  =  5.5 

0.0000] [10.263 1  26.2500  0.0000][5.5000  29.0000 

r\  =  (l/10)tw=10/10=l 

0.0000][0.0000  23.5000  0.0000] 

a  =  nt  6 

hf  =  22.82 

Af  =  490.54 

Table  27  Data  for  Flange  Model  8  in  Table  2  (units  in  mm) 


6,  =  7 

[0.0000  0.7000  0.0000][0.2050  0.2050  0.0000][0.7000  0.0000 

f  =22 

0.0000][25.6299  0.0000  0.0000][25.6299  7.0000 

1 bf 

0.0000] [22.6299  7.8038  0.0000][20.4338  10.0000 

r  —  6 

0.0000][1 1.1962  26.0000  0.0000][6.0000  29.0000 

II 

o 

s 

II 

o 

0.0000][0.0000  23.0000  0.0000] 

a-n  /  6 

hf  =25.63 

4  =  519.51 
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Table  28  Data  for  Flange  Model  10  in  Table  2  (units  in  mm) 


K=  9 

[0.0000  0.9000  0.0000][0.2636  0.2636  0.0000][0.9000  0.0000 

t  =22 

0.0000][25.6299  0.0000  0.0000][25.6299  9.0000 

lbf 

0.0000] [22.6299  9.8038  0.0000][20.4338  12.0000 

r  —  6 

0.0000][1 1.1962  28.0000  0.0000][6.0000  31.0000 

II 

o 

*-~+ 

s 

II 

vO 

O 

0.0000][0.0000  25.0000  0.0000] 

a  =  nt  6 

hf  =25.63 

4  =  570.70 

Table  29  Data  for  Flange  Model  12  in  Table  2  (units  in  mm) 


t  =8 

w 

[0.0000  0.8000  0.0000][0.2343  0.2343  0.0000][0.8000  0.0000 

/  =  25 

0.0000][29.5167  0.0000  0.0000][29.5167  8.0000 

bf 

0.0000][26.0167  8.9378  0.0000][23.4545  11.5000 

r  =  7 

0.0000][13.0622  29.5000  0.0000][7.0000  33.0000 

o 

00 

II 

* 

o' 

II 

0.0000][0.0000  26.0000  0.0000] 

a  =  ni  6 

hf  =  29.52 

Af  =  681.94 

Table  30  Data  for  Flange  Model  14  in  Table  2  (units  in  mm) 


K  =  10 

[0.0000  1.0000  0.0000][0.2929  0.2929  0.0000][1.0000  0.0000 

/  =  25 

0.0000][29.5167  0.0000  0.0000][29.5167  10.0000 

0.0000][26.0167  10.9378  0.0000][23.4545  13.5000 

r  =  7 

0.0000][13.0622  31.5000  0.0000][7.0000  35.0000 

r\  =  (l/10)tw=10/10=l 

0.0000][0.0000  28.0000  0.0000] 

a  =  ;r/6 

=  29.52 

Af  =  740.89 
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2. 


DATA  FOR  THE  BULB-FLAT  STIFFENER  CROSS-SECTION 


See  Figure  1,  Figure  7,  and  Figure  8. 


Table  31  Data  for  first  stiffener  in  Table  5  (units  in  mm) 


o 

<N 

II 

-s:' 

[0.0000  0.0000  0.0000][6.0000  0.0000  0.0000][6.0000  99.4115 

hr  =20  59 

0.0000][6.6699  101.9115  0.0000][8.5000  103.7417 

0.0000][20.5000  1 10.6699  0.0000][23.0000  115.0000 

K  =99.41 

0.0000] [1 8.0000  120.0000  0.0000][0.6000  120.0000 

K  =  6 

0.0000] [0. 1757  119.8243  0.0000][0.0000  119.4000  0.0000] 

V=17 

r  =  5 

II 

o 

s 

II 

ON 

o 

a  =  ?r/6 

Table  32  Data  for  second  stiffener  in  Table  5  (units  in  mm) 


hs  =  120 

[0.0000  0.0000  0.0000][8.0000  0.0000  0.0000][8.0000  99.4115 

hr  =20  59 

0.0000][8.6699  101.9115  0.0000][10.5000  103.7417 

0.0000][22.5000  1 10.6699  0.0000][25.0000  115.0000 

hw  =99.41 

0.0000] [20.0000  120.0000  0.0000][0.8000  120.0000 

C=8 

0.0000][0.2343  1 19.7657  0.0000][0.0000  119.2000  0.0000] 

v=17 

r  =  5 

o 

00 

II 

£ 

o" 

II 

5C" 

a  =  n !  6 
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Table  33  Data  for  third  stiffener  in  Table  5  (units  in  mm) 


hs  =  140 

[0.0000  0.0000  0.0000][7.0000  0.0000  0.0000][7.0000  117.1795 

h  =  22.82 

0.0000][7.7369  119.9295  0.0000][9.7500  121.9426 

0.0000][23.2500  129.7369  0.0000][26.0000  134.5000 

K  =117.18 

0.0000][20.5000  140.0000  0.0000][0.7000  140.0000 

K  =  7 

0.0000][0.2050  139.7950  0.0000][0.0000  139.3000  0.0000] 

v=19 

r  =  5.5 

II 

o 

s 

II 

o 

a-n  /  6 

Table  34  Data  for  fourth  stiffener  in  Table  5  (units  in  mm) 


hs  =  140 

[0.0000  0.0000  0.0000][10.0000  0.0000  0.0000][10.0000 

h  =  22  82 

117.1795  0.0000][10.7369  119.9295  0.0000][12.7500  121.9426 

0.0000][26.2500  129.7369  0.0000][29.0000  134.5000 

K  =117.18 

0.0000][23.5000  140.0000  0.0000][1.0000  140.0000 

6,  =  10 

0.0000] [0.2929  139.7071  0.0000][0.0000  139.0000  0.0000] 

V=19 

r  =  5.5 

6  =  (l/10)tw=10/10=l 

a-n  /  6 

Table  35  Data  for  fifth  stiffener  inTable  5  (units  in  mm) 


hs  =  160 

[0.0000  0.0000  0.0000][7.0000  0.0000  0.0000][7.0000  134.3701 

hr  =25  63 

0.0000][7.8038  137.3701  0.0000][10.0000  139.5662 

0.0000] [26.0000  148.8038  0.0000][29.0000  154.0000 

K  =134.37 

0.0000][23.0000  160.0000  0.0000][0.7000  160.0000 

tw  =  7 

0.0000][0.2050  159.7950  0.0000][0.0000  159.3000  0.0000] 

tbj  =  22 

r  =  6 

ii 

o 

s 

II 

o 

a-n /6 
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Table  36  Data  for  sixth  stiffener  in  Table  5  (units  in  mm) 


hs  =  160 

[0.0000  0.0000  0.0000][9.0000  0.0000  0.0000][9.0000  134.3701 

h,  =  25  63 

0.0000]  [9.8038  137.3701  0.0000][12.0000  139.5662 

0.0000][28.0000  148.8038  0.0000] [3 1.0000  154.0000 

K  =134.37 

0.0000][25.0000  160.0000  0.0000][0.9000  160.0000 

K=  9 

0.0000][0.2636  159.7364  0.0000][0.0000  159.1000  0.0000] 

<N 

<N 

II 

r  =  6 

II 

o 

s 

II 

vO 

O 

a  =  n  /  6 

Table  37  Data  for  seventh  stiffener  in  Table  5  (units  in  mm) 


o 

00 

II 

-sr 

[0.0000  0.0000  0.0000][8.0000  0.0000  0.0000][8.0000  150.4833 

h  =  29.52 

0.0000] [8.9378  153.9833  0.0000] [  1 1 .5000  156.5455 

0.0000][29.5000  166.9378  0.0000][33.0000  173.0000 

hw  =150.48 

0.0000] [26.0000  180.0000  0.0000][0.8000  180.0000 

*w=8 

0.0000][0.2343  179.7657  0.0000][0.0000  179.2000  0.0000] 

<N 

II 

r  =  1 

jyt 

II 

o 

*-+. 

II 

00 

o 

a  =  n  /  6 

Table  38  Data  for  eighth  stiffener  in  Table  5  (units  in  mm) 


o 

00 

II 

[0.0000  0.0000  0.0000][10.0000  0.0000  0.0000][10.0000 

hf  =  29.52 

150.4833  0.0000][10.9378  153.9833  0.0000][13.5000  156.5455 

0.0000][31.5000  166.9378  0.0000][35.0000  173.0000 

hw  =150.48 

0.0000][28.0000  180.0000  0.0000][1.0000  180.0000 

C  =  io 

0.0000] [0.2929  179.7071  0.0000][0.0000  179.0000  0.0000] 

hf  =  25 

r  =  1 

^=(1/10)^=10/10=1 

a  =  n  /  6 
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Table  39  Data  for  angle  stiffener  in  Figure  9  (units  in  mm) 


K  =  120 

[0  0  0] [6  0  0] [6  107.5259068  0][23  107.5259068  0][23  120  0][0  120  0] 

Kva  ~  107.53 

K  =  6 

v=17 

<n 

II 

3.  DATA  FOR  THE  STIFFENED  PLATE  MODELS 


Table  40  Data  for  stiffened  plate  models  1-3  with  no  flange  (units  in  inches) 


See  Table  7  for 
parameters. 

[0.0000  0.0000  3.3233][72.0000  0.0000  3.3233][0.0000  0.0000 
0.0000] [72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
-20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 

Table  41  Data  for  stiffened  plate  models  4-6  with  no  flange  (units  in  inches) 

See  Table  10  for 
parameters. 

[0.0000  0.0000  3.9138][72.0000  0.0000  3.9138][0.0000  0.0000 
0.0000] [72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
-20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 

Use  the  following  tables  to  build  the  stiffened  plate  models.  In  MSC  Patran, 
create  the  geometry  from  the  points  in  each  table  by  entering  the  points  and  then  create 
the  curves  and  surfaces  as  appropriate.  Create  a  mesh  (with  mesh  size  1)  for  the  plate  and 
web  surfaces.  Create  a  mesh  for  the  curve  at  the  top  of  the  web  with  topology  BAR2. 
The  curve  will  serve  as  the  flange  once  the  appropriate  beam  section  is  created.  Using 
the  properties  menu,  create  and  name  a  ID  beam  property  set.  Use  the  general  section 
(CBEAM)  option.  Use  the  input  properties  button  and  then  the  create  sections  beam 
library  to  create  and  emplace  the  appropriate  flange  onto  the  web.  For  the  bulb-flat 
flange,  use  the  beam  library  to  create  and  name  a  bulb-flat  flange  by  selecting  the 
arbitrary  shape  button  and  the  boundary  loops  method  button.  For  the  input  option,  set 
the  select  surface  option.  In  the  select  surface  menu,  set  the  maximum  allowable 
curvature  error  to  0.005.  Select  the  surface  and  press  OK.  In  the  general  beam 
(CBEAM)  input  properties  subordinate  menu,  enter  the  appropriate  material  name, 
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section  name  (just  created  in  the  create  sections  beam  library),  and  bar  orientation  (start 
by  using  the  vector  (0, 1, 0)  and  adjust  with  (l,0,0)  or  (0, 0,  l)  as  needed). 


Table  42  Data  for  bulb-flat  stiffened  plate  model  lb  (units  in  inches) 


See  Table  7  model 
lb  for  parameters. 

[0.0000  -0.1181  3.3233][0.0000  0.1181  3.3233][0.0000  0.1445 
3.4217][0.0000  0.2165  3.4938][0.0000  0.6890  3.7665][0.0000 

0.7874  3.9370][0.0000  0.5906  4.1339][0.0000  -0.0945 
4.1339][0.0000  -0.1 1 12  4.1269][0.0000  -0.1 181  4.1 102][0.0000 
0.0000  3.3233][72.0000  0.0000  3.3233][0.0000  0.0000 

0.0000] [72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
-20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 

Table  43  Data  for  bulb-flat  stiffened  plate  model  2b  (units  in  inches) 

See  Table  8  model 
2b  for  parameters. 

[0.0000  -0.1181  3.3233][0.0000  0.1181  3.3233][0.0000  0.1445 
3.4217][0.0000  0.2165  3.4938][0.0000  0.8071  3.8347][0.0000 

0.9055  4.0052][0.0000  0.7087  4.2020][0.0000  -0.0945 

4.2020] [0.0000  -0.1 1 12  4. 1 95 1  ] [0.0000  -0.1181  4.1784][0.0000 
0.0000  3.3233][72.0000  0.0000  3.3233][0.0000  0.0000 

0.0000] [72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
-20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 

Table  44  Data  for  bulb-flat  stiffened  plate  model  3b  (units  in  inches) 

See  Table  9  model 
3b  for  parameters. 

[0.0000  -0.1181  3.3233][0.0000  0.1181  3.3233][0.0000  0.1445 
3.4217][0.0000  0.2165  3.4938][0.0000  1.0039  3.9484][0.0000 

1.1024  4.1 189][0. 0000  0.9055  4.3157][0.0000  -0.0945 

4.3 157][0. 0000  -0.1 1 12  4.3088][0.0000  -0.1181  4.2921][0.0000 
0.0000  3.3233][72.0000  0.0000  3.3233][0.0000  0.0000 

0.0000] [72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
-20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 

Table  45  Data  for  bulb-flat  stiffened  plate  model  4b  (units  in  inches) 

See  Table  10 
model  4b  for 
parameters. 

[0.0000  -0.1181  3.9138][0.0000  0.1181  3.9138][0.0000  0.1445 
4.0123][0.0000  0.2165  4.0843][0.0000  0.6890  4.3571][0.0000 

0.7874  4.5276][0.0000  0.5906  4.7244][0.0000  -0.0945 

4.7244] [0.0000  -0.1 1 12  4.7175][0.0000  -0.1 181  4.7008][0.0000 
0.0000  3.9138][72.0000  0.0000  3.9138][0.0000  0.0000 

0.0000] [72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
-20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 
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Table  46  Data  for  bulb-flat  stiffened  plate  model  5b  (units  in  inches) 

See  Table  1 1  [0.0000  -0. 1  1 8 1  3.9 1 38][0.0000  0. 1  1 8 1  3.9 1 38][0.0000  0. 1 445 

model  5b  for  4.0 123] [0.0000  0.2165  4.0843][0.0000  0.8071  4.4253][0.0000 

parameters.  0.9055  4.5958] [0.0000  0.7087  4.7926][0.0000  -0.0945 

4.7926] [0.0000  -0.1 1 12  4.7857][0.0000  -0.1 181  4.7690][0.0000 
0.0000  3.9138][72.0000  0.0000  3.9138][0.0000  0.0000 
0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
_ -20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 


Table  47  Data  for  bulb-flat  stiffened  plate  model  6b  (units  in  inches) 

See  Table  12  [0.0000  -0. 1  1 8 1  3.9 1 38][0.0000  0. 1  1 8 1  3.9 1 38][0.0000  0. 1 445 

model  6b  for  4.0123][0.0000  0.2165  4.0843][0.0000  1.0039  4.5389][0.0000 

parameters.  1.1024  4.7094] [0.0000  0.9055  4.9063][0.0000  -0.0945 

4.9063][0.0000  -0.1 1 12  4.8993][0.0000  -0.1 181  4.8826][0.0000 
0.0000  3.9138][72.0000  0.0000  3.9138][0.0000  0.0000 
0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
_ -20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 


To  build  the  circular  flange  stiffened  plates,  ensure  a  mesh  is  created  for  the  curve 
at  the  top  of  the  web  with  topology  BAR2.  The  curve  will  serve  as  the  flange  once  the 
appropriate  beam  section  is  created.  Using  the  properties  menu,  create  and  name  a  ID 
beam  property  set.  Use  the  general  section  (CBEAM)  option.  Use  the  input  properties 
button  and  then  the  create  sections  beam  library  button  to  create  and  emplace  the 
appropriate  flange  onto  the  web.  Use  the  beam  library  to  create  and  name  a  flange  by 
selecting  the  standard  shape  cross-section  using  the  solid  circular  shape  button  with  the 
appropriate  radius  dimension.  In  the  general  beam  (CBEAM)  input  properties 
subordinate  menu,  enter  the  appropriate  material  name,  section  name  (just  created  in  the 
create  sections  beam  library),  and  bar  orientation. 

Table  48  Data  for  circular  flange  stiffened  plate  model  lc  (units  in  inches) 

See  Table  7  model  [0.0000  0.0000  3.730 1 97][72.0000  0.0000  3.730 1 97][0.0000 
lc  for  parameters.  0.0000  0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000 
0.0000][72.0000  -20.0000  0.0000][0.0000  20.0000 
_ 0.0000] [72.0000  20.0000  0.0000] _ 
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Table  49  Data  for  circular  flange  stiffened  plate  model  2c  (units  in  inches) 


See  Table  8  model 
2c  for  parameters. 

[0.0000  0.0000  3.768468][72.0000  0.0000  3.768468][0.0000 

0.0000  0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000 

0.0000] [72.0000  -20.0000  0.0000][0.0000  20.0000 

0.0000] [72.0000  20.0000  0.0000] 

Table  50  Data  for  circular  flange  stiffened  plate  model  3c  (units  in  inches) 

See  Table  9  model 
3c  for  parameters. 

[0.0000  0.0000  3.831455][72.0000  0.0000  3.831455][0.0000 

0.0000  0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000 

0.0000] [72.0000  -20.0000  0.0000][0.0000  20.0000 

0.0000] [72.0000  20.0000  0.0000] 

Table  51  Data  for  circular  flange  stiffened  plate  model  4c  (units  in  inches) 

See  Table  10 
model  4c  for 
parameters. 

[0.0000  0.0000  4.320748][72.0000  0.0000  4.320748][0.0000 

0.0000  0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000 

0.0000] [72.0000  -20.0000  0.0000][0.0000  20.0000 

0.0000] [72.0000  20.0000  0.0000] 

Table  52  Data  for  circular  flange  stiffened  plate  model  5c  (units  in  inches) 

See  Table  1 1 
model  5  c  for 
parameters. 

[0.0000  0.0000  4.359019][72.0000  0.0000  4.359019][0.0000 

0.0000  0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000 

0.0000] [72.0000  -20.0000  0.0000][0.0000  20.0000 

0.0000] [72.0000  20.0000  0.0000] 

Table  53  Data  for  circular  flange  stiffened  plate  model  6c  (units  in  inches) 

See  Table  12 
model  6c  for 
parameters. 

[0.0000  0.0000  4.422005][72.0000  0.0000  4.422005][0.0000 

0.0000  0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000 

0.0000] [72.0000  -20.0000  0.0000][0.0000  20.0000 

0.0000] [72.0000  20.0000  0.0000] 

To  build  the  T-flange  stiffened  plate  models,  the  same  data  in  Table  54  is  used  to 
build  T-flange  stiffened  plate  models  It,  2t,  and  3t.  The  same  data  in  Table  55  is  used  to 
build  T-flange  stiffened  plate  models  4t,  5t,  and  6t.  Ensure  a  mesh  is  created  for  the 
curve  at  the  top  of  the  web  with  topology  BAR2.  The  curve  will  serve  as  the  flange  once 
the  appropriate  beam  section  is  created.  Using  the  properties  menu,  create  and  name  a 
ID  beam  property  set.  Use  the  general  section  (CBEAM)  option.  Use  the  input 
properties  button  and  then  the  create  sections  beam  library  button  to  create  and  emplace 
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the  appropriate  flange  onto  the  web.  In  the  beam  library,  create  and  name  a  standard 
shape  cross-section  using  the  solid  rectangular  shape  button  with  the  appropriate 
dimensions.  In  the  general  beam  (CBEAM)  input  properties  subordinate  menu,  enter  the 
appropriate  material  name,  section  name  (just  created  in  the  create  sections  beam  library), 
and  bar  orientation  (start  by  using  the  vector  (0,1,0)  and  adjust  with  (l,0,0)  or  (0,0,l) 
as  needed). 


Table  54  Data  for  T-flange  stiffened  plate  models  It,  2t,  and  3t  (units  in  inches) 


See  Table  7  model 

It,  Table  8  model 

2t,  or  Table  9 
model  3t  as 
appropriate  for 
parameters. 

[0.0000  0.0000  3.4414][72.0000  0.0000  3.4414][0.0000  0.0000 
0.0000] [72.0000  0.0000  0.0000][0.0000  -20.0000  0.0000][72.0000 
-20.0000  0.0000][0.0000  20.0000  0.0000][72.0000  20.0000  0.0000] 

Table  55  Data  for  T-flange  stiffened  plate  model  4t,  5t,  and  6t  (units  in  inches) 

See  Table  10 
model  4t,  Table  1 1 
model  5t,  or  Table 

12  model  6t  as 
appropriate  for 
parameters. 

[0.0000  0.0000  4.03 195 1][72.0000  0.0000  4.03 195 1][0.0000 

0.0000  0.0000][72.0000  0.0000  0.0000][0.0000  -20.0000 

0.0000] [72.0000  -20.0000  0.0000][0.0000  20.0000 

0.0000] [72.0000  20.0000  0.0000] 
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